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ABSTRACT 


The reaction concept is employed to formulate an integral 
equation for radiation and scattering from plates, corner reflectors, 
and dielectric-coated conducting cylinders. The surface-current 
density on the conducting surface is expanded with subsectional bases. 
The dielectric layer is modeled with polarization currents radiating 
in free space. Maxwell’s equation and the boundary conditions are 
employed to express the polarization-current distribution in terms of 
the surface-current density on the conducting surface. By enforcing 
reaction tests with an array of electric test sources, the moment 
method is employed to reduce the integral equation to a matrix equation. 
Inversion of the matrix equation yields the current distribution, and 
the scattered field is then obtained by integrating the current 
distribution. 

This report presents the theory, computer program and numerical 
results for radiation and scattering from plates, corner reflectors, 
and dielectric-coated conducting cylinders. 
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CHAPTER I 
INTRODUCTION 


Low-frequency solutions for radiation and scattering from 
cylinders have been reported in several published papers [l-8] . Among 
them, Mei and Van Bladel |_l] employed a point-matching procedure to 
solve the electric-field integral equation and the magnetic-field 
integral equation for transverse-magnetic and transverse-electric 
incident waves, respectively. In the point-matching procedure the 
surface current distribution was expanded into rectangular-pul se bases 
and the appropriate boundary conditions were enforced at discrete 
points on the conducting surface of the cylinder. Richmond [4] has 
developed a wire-grid array model for cylinders with transverse- 
magnetic incident wave and shown that if a sufficiently great number of 
wires is employed, the scattering pattern approaches that of a solid 
cylinder of the same contour. Richmond [9] has also developed a 
piecewise-sinusoidal reaction formulation for electromagnetic radiation 
and scattering problems involving cylinders with non-circular cross 
section for the transverse-electric incident wave. No solution, 
however, has been published for non -circular cylinders with a 
dielectric coating. 

For three-dimensional problems, two methods are available for 
electromagnetic modeling of a continuous conducting surface with 
arbitrary shape: the wire-grid model [10] and the surface-current 

model [11,12] with rectangular-pulse bases. Both methods have similar 
limitations with the maximum cell width restricted to approximately 
one-tenth of a wavelength. Unless the conducting body is symmetric or 
is a figure of revolution, computer storage requirements have limited 
the conducting surface area to one or two square wavelengths. 

In this dissertation, Rumsey's [13] reaction concept is employed 
to formulate an integral equation for scattering by plates, comer 
reflectors and dielectric-coated cylinders with noncircular cross 
section. In the reaction formulation, the surface -current density 
on the conducting surface is expanded with suitable bases. The 
dielectric layer is modeled with the equivalent polarization currents 
radiating in free space. Maxwell's equations and the boundary con- 
ditions are employed to express the polarization-current distribution 
in terms of the surface-current density on the conducting surface. 

By enforcing reaction tests with an array of electric test sources, 
the moment method is employed to reduce the integral equation to a 
matrix equation. Numerical solution of this system yields a station- 
ary result for the samples of the current distribution. Finally, 
the quantities of interest such as the gain, far-field pattern, and 
the radar cross section are determined from the current distribution. 
With perfect conductivity, the analysis presented in this dissertation 
is valid for open as well as closed cylinders. With finite conduc- 
tivity or with a thin dielectric coating, however, the analysis is 
restricted to closed cylinders. 
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The remaining text presents the general theory of the reaction 
formulation for radiation and scattering from conducting bodies. The 
time dependence ed^^ is understood and suppressed. Chapter II pre“ 
sents the detailed theoretical outline of the reaction concept which 
forms the foundation of this dissertation. Relevant electric sources 
employed as the expansion functions and the test sources are discussed 
in Chapter III. Evaluation of the mutual impedances and excitation 
voltages are considered in Chapter IV and Chapter V. The electro- 
magnetic modeling of the conducting body is an essential step in the 
reaction formulation. This is described in Chapter VI. The field 
scattered by a conducting body is the sum of the contributions from 
all the current modes which appear in the expansion for the surface- 
current distribution. In Chapter VII the far-fi eld contributions of 
these sources are discussed. Numerical results for the radar cross 
section of rectangular plates, comer reflectors and dielectric- 
coated cylinders are presented in Chapter VIII. 
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CHAPTER II 

THE REACTION INTEGRAL FORMULATION 


Two well known integral equations, the electric-field integral 
equation and the magnetic-field integral equation are usually employed 
to solve electromagnetic radiation and scattering problems. In this 
chapter, however, we develop the more general reaction integral equa- 
tion of Rumsey [l3j . It has been noted [14] that the reaction inte- 
gral equation is more general in the sense that it can be reduced 
to either the electric-field integral equation or the magnetic- field 
integral equation if one enforces the reaction integral equation with 
a set of delta-function electric or magnetic test sources. In the 
following sections the reaction concept and its application in 
electromagnetic problems will be examined. 

Consider the exterior scattering problem illustrated in Fig. la. 
In the presence of a conducting body, the impressed electric and mag- 
netic currents (Oi ,Mi ) generate the electric and magnetic field 
intensities (£,Hj. For simplicity, let the exterior medium be free 
space. 

From the surface-equivalence theorem of Schelkunoff [l5j , the 
interior field will vanish (without disturbing the exterior field) if 
we introduce the following surface-current densities 


(1) Jg = n X H^ 

(2) Ms " L ^ ^ 

on the closed surface S of the scatterer. (The unit vector n is 
directed outward on S.) In this situation, illustrated in Fig. lb, we 
may replace the scatterer with free space without disturbing the 
field anywhere. 

By definition, the incident field (EjiHi) generated by 
in free space, and the scattered field is: 


( 3 ) h = 


(4) Hg = H - Hi . 


When the surface current (^>M^) radiates in free space, it generates 
the field (^,Hc) in the e^enor and (-Ej,-H^-) in the interior region 
Tbis result, illustrated in Fig. Ic, is deduced from Fig. lb and the 
superposition theorem. 
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Fig. la— The source generates the field (^»^) 

with scatterer. 



Fig. lb — The interior field vanishes when the currents 
are introduced on the surface of the 
scatterer. 


r" 

'Vv FREE SPACE ' 

K y 




Fig. Ic — The exterior scattered field may be generated 
by in free space. 
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With the scatte rer replaced by free space, we have noted in 
Fig. lb that the interior region has a null field. As shown in Fig. 2, 
we place an electric test source ^ in this region and find from the 
reciprocity theorem that 

(5) 1^5 (i • it - Ms • Ht) ds + fff (J^ • ^ - Mi . Ht) dv = 0 

where is the free-space field of the test source. In words, 

Eq. (5) stares that the interior test source has zero reaction with the 
other sources. This "zero- re action theorem" was developed by 
Rumsey \ l3] . 

Equation (5) is the integral equation for the scattering problem, 
and our objective is to use this equation to determine the surface - 
current distributions ^ and To accomplish this, we expand these 
functions in finite senes so there will be a finite number N of un- 
known expansion constants. Next we obtain N simultaneous linear 
equations to permit a solution for these constants. One such equation 
is obtained from Eq. (5) each time we set up a new test source. 





V^\FREE SPACE 



Fig. 2— An electric test source is positioned in 
the interior of the scattering region. 


The magnetic current vanishes if the scatterer is a perfect 
conductor. We assume a finite conductivity and use the impedance 
boundary condition: 


(6) = Zg 4 x n 

where Zg denotes the surface impedance defined by n x £ = Zg n x 
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For three-dimensional problems involving arbitrary scatter^rs, 
^ and are functions only of the position on the surface of the 
scatterer. Eqs. (5) and (6) yield 


C7) 


-f 


JJ 






- (n X |i„) Zsl ds = Jj J. 


Em 


ds 





where (E^^>H^) denotes the free-space field of test-source m. 

We represent the electric current distribution as follows: 


N 

(8) i = Z In i 

n=l 


where the complex constants 1,^ are samples of the function J The 
vector functions are known as basis functions, subsectional bases, 
expansion functions or dipole modes. We employ expansion functions 
^ and test sources ^ with unit current density at the terminals. 

From Eqs. (7) and (8) we obtain the simultaneous linear equations 

N 

Z In \ = 1>2,3, N 

n=l 

where 


(10) = -jj^ ^ • [E^ - (n X ds — Jj ^ • E^ ds 

01) = j] ds -jj . Mi . ds = JJ . E. ds 

In Eqs. (10) and (11) the integrations extend over the region where 
the integrand is non-zero. For example, region n is that portion 
of the surface S covered by the expansion function Region m 
covers the interior test source The reciprocity theorem relates 
the first and second integrals in Eq. (10). In the second integral, 
^ is the free-space field generated by ^ and the associated 
magnetic current 

For computational speed and storage, it will be advantageous 
to have a symmetric impedance matrix Furthermore, the test 

sources should be selected to yield a wel 1-condi tioned set of 
simultaneous linear equations. For these reasons, we employ test 
sources ^ of the same size, shape and functional form as the 
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expansion functions Final we position the interior test sources 
a small distance 6 from surface S and take the limiting form of the 
integrals as s tends to zero. 

The effect of a dielectric coating on a conducting body will 
now be considered. For simplicity, let the dielectric layer have 
the same permeability as free space. From the volume equivalence 
theorem of Rhodes [16] the dielectric coating may be replaced with 
free space and an equivalent electric current density 


(12) = j4e-eo) E 


where E denotes the electric field intensity in the dielectric and 
e = is the permittivity of the dielectric layer. From Eq. (12) 
the equivalent current J vanishes outside the region of the 
dielectric coating. ^ 

Let (£,H^) denote the field generated by in the presence 

of a dielectric-coated conducting scatterer. (Tutside the scatterer, 
this field may also be generated by and radi- 

ating in free space. These sources, radFating in free space, 
generate a null field in the interior region of the conducting body. 
The surface currents (^,^) are located on the surface of the 
conducting body and are related to the field (E_,h[) by Eqs. (1) and 


For a coated conducting body, the reaction integral equation 
(Eq. (7)) is modified by replacing with + ^q. The current 
density may be regarded as an additional source which plays 
much the same role as the impressed source J-j . However, is an 
unknown quantity because ^ is unknown. If the dielectric coating is 
thin. Maxwell's equations and boundary conditions can be employed 
to express the polarization-current distribution in terms of the 
surface -current density on the conducting surface. ^Therefore, 
may be regarded as a dependent unknown function because it is 
simply related to 

The polarization- current distribution is expanded as follows: 

iq ~ I ^n iin 
n=l 

where the ^ are functions simply related to Thus, for a 
coated conducting body, each expansion mode in Eq. (8) has 
associated with it a polarization current and the reaction 
between the electric test source m and the expansion mode n has an 
additional term given by 
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( 14 ) 


AC 


mn 




Ijn 


ds 


where the integration extends through the dielectric coating in the 
range of the expansion mode n. The functions ^ are defined over a 
surface, while the ^ are defined in a volumetric region. 

It may be noted that in the reaction formulation, the effects 
of a dielectric coating are accounted for entirely through a modi- 
fication of the square reaction matrix This modification 

influences the current distribution, field patterns and scattering 
properties. 

The following chapter discusses the electric sources which ar^ 
employed as test sources and expansion modes for the current 
distribution on the conducting surface. 
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CHAPTER III 

ELECTRIC TEST SOURCES AND EXPANSION MODES 


In the preceding chapter, the reaction integral equation is 
reduced to a matrix equation via the moment method in two steps. 

First the unknown current distribution is expanded with basis func- 
tions, Test sources with the same functional form are then employed 
to perform the reaction tests. In the following sections three types 
of test sources employed in this dissertation are discussed. 

A, LONGITUDINAL STRIP SOURCES 

For an infinitely long cylinder with transverse magnetic incident 
wave, the electric current density induced on the conducting surface 
is in the longitudinal direction. Thus a suitable choice of the 
bases is the longitudinal strip source. 

Consider the "strip source" illustrated in Fig. 3. This source 
is an electric surface -current distribution J = z J(x) located on the 
xz plane. This source has width h and infinite length and radiates 
in free space. For the strip source shown in Fig. 3, the fields are 
given as follows: 


A 



Fig. 3— An electric strip source and the coordinate 
system. 


(15) 



0 Hj2)(kp) dx' 
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(16) 


- iiE 1 
4 j 

^ X p Hj^\kp) 

0 

dx' 

where 





(17) 

p = 

/{x-x* 

-yrvY^ 


(18) 

k = 




(19) 

n = 




If 

the 

electric strip source 

has a uniform distribution. 

(20) 

J(x) 

A 

= z 




Eqs. (15) and (16) can be written as 


(21) 

r _ 

^ ■ 4 

<"2) , 

Hq (kp) dx' 

^0 

(22) 

J 

z X p (kp) 

0 


B. TRANSVERSE STRIP SOURCES [9] 

In the case of a conducting cylinder illuminated by a transverse- 
electric incident wave, the electric current induced on the surface 
is in the transverse direction. Therefore, transverse electric 
sources are the natural choice for the induced current density. 

Again consider the’^trip source" illustrated in Fig. 3. But this 
time the source is an electric surface-current distribution ^ = x J(x) 
located on the xz plane. For this transverse electric source the 
free space fields are 

(23) E = -;r (kp) dx' +T 2 J' h' (kp) dx' 

■'0 ^ 0 

(24) H = - I J X ? Hi^’(kp) dx- 
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where 


( 25 ) J' = dJ/dx' 


For most current functions J(x), the field integrals must be 
evaluated with infinite-series expansions or numerical integration 
procedures. For the sinusoidal current distribution, however, Ew 
is obtained in simple closed form [ 9 ]. Thus, if 

( 26 ) " x|IjSin(kh-kx) + l2sin(kx)]/sin(kh) 


then 

DiHo"’(kpi) cos(kh) - I,Ho“Vp2)J 

^ Sin(kh) 

+ [l2Hp^\kp2) cos(kh) - I2 Hq^^Upi)J 

where Ij and K represent 0(0) and 0(h), respectively, and pn, p^ 
are illustrated in Fig. 3 . 

For the purpose of representing a continuous surface -current 
distribution on a conducting cylinder, it is useful to define a 
strip dipole which is corpprised of two strip monopoles. Fig. 4 a 
illustrates a planar strip dipole. This dipole lies on the xz plane 
and has infinite length in the z direction. The surface -current 
density is 

^ sin k(x-Xi) 

^ X, < x< X2 

sin k(x2“Xj) 


^ sin k()<3 -x) 

( 29 ) “ X for X2< x< x 

sin k()^-X2) ^ 


As indicated in Fig. 4 b, the current density vanishes at the edges 
Xi and X3 , is continuous across the terminals at X2 and has a slope 
discontinuity at X2. 
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y 


J j 



X| Xg X 3 

Fig. 4b — The current-densi ty distribution ^ on the 
sinusoidal strip dipole. 


Fig. 5 illustrates a strip V-dipole. Distance along the dipole 
arms is measured by the coordinates s and t with origin at the 
terminals 0. The surface -current density is 

^ sin k(si-s) 

(30) ^ =-s — — on arm s 

sin ksj 

^ sin k(tj-t) 

(31) i = t — — — - — on arm t, 

sin ktj 


where the unit vectors s and t are perpendicular to the z-axis. Thus, 
the current density vanishes at edges s. and t^ and has unit value 
at the terminals 0. The edges are parallel with the z-axis. If the 
wedge angle is adjusted to 180 degrees, the V-dipole in Fig. 5 
reduces to the planar dipole in Fig. 4. 
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Fig. 5“-Nonplanar strip dipole with edges at Sj and 
tj and terminals at 0. 


C. RECTANGULAR SURFACE SOURCES 


Consider a "surface monopole " radiating in free space as shown 
in Fig. 6. This source is an electric surface -current density 
^ = z J(z) located on the yz plane. This source has height c(/2 and 
width b. The surface -current density is related to the current by 
I^= b J. For the electric surface monopole illustrated in Fig. 6, 
the fi^ds are 


(32) 



r^/2 ,b 

K2 

L •'o JO 


J G^(kR) dy' 


dz' +v 


a/2 

0 


b "1 

J'GjkR)dy'dz' 
Jo 


,a/2 b 

(33) H = V x[ ( J G (kR) dy' dz‘ 
Jo Jo ° 


where 


(34) R = /Tx^^^“(y-y'y + (z-z')^ 

(35) J' = dJ/dz* 
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(x.y.z) 



Fig. 6- -An electric surface inonopole and the coordinate 
system. 


(36) Go(kR) = 


Numerical integration techniques must be employed to perform the 
fiel d integrals . 


A planar sinusoidal dipole source located on the yz plane as 
shown in Fig. 7a will be considered. This source is an electric 
surface-current density with height a and width b. The surface- 
current density is given by 


(37) 

(38) 


0 = z cos 


/ (z-zj- 


V2(Z2~Zi) 


0 = 


z cos 


' (z-Z2)tt'\ 

. 2(23 


for Zj< z< Z2 


for z^^ Z5i Z3 


As illustrated in Fig. 7b and 7c, the current density vanishes at 
the edges z = Zj and z = Z3 , and is uniformly distributed in the 
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2 Z 



Fig. 7— An electric surface dipole and its current- 
density distribution. 


transverse direction. The surface -current density and its slope are 
continuous across the terminals at z - Z 2 for ~ a/2. 

Fig. 8 illustrates a surface V-dipole. Distance along the 
dipole arms is measured by the coordinates s and t with origin at the 
terminal 0. The surface -current density for this electric source is 

(39) £ = - s cosf J on arm s 
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on arm t . 


(40) J = t cos 


When the wedge angle 4 / is adjusted to 180 degrees, the V-dipole 
in Fig , 8 reduces to the planar surface dipole in Fig. 7a. 


A 

t 

/ 



A 

S 



Fig. 8 --A nonplanar surface dipole with edges at s, and 
tj and terminal at 0 . 


The electric sources defined in the previous sections are all 
hypothetical sources. The current density on a conducting strip is 
neither uniform nor sinusoidal and the current density induced on a 
rectangular plate is not sinusoidal. The relevance of these sources 
will now be explained. The electric current distributions (Eqs. ( 20 ), 
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(30), (31), (39), and (40)) will be used as the basis functions 
(Eq. (8)) for expanding the unknown current distribution induced on 
conducting surfaces for various problems. Furthermore, test sources 
with the same size, shape, and functional fonn as the expansion 
functions will be employed with the reaction concept to solve the 
integral equation. 

By superposition, the field scattered by a conducting body may 
be regarded as the sum of the fields radiated by the mode currents: 

^ n=i ^ 

W h" = I In ^ 

n=l 

where (e^^, h^) is the free space field generated by the mode current 

n , 
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CHAPTER IV 
■mE IMPEDANCE MATRIX 


From the viewpoint of reaction, the complex number C in 
Eq. (10) represents the reaction between the sources m anTn. The 
)^action between two electric sources m and n is: 


(43) 


'mn 



in 


Although the electric sources defined in Chapter III are 
hypothetical, it is useful to define self impedance with the induced- 
emf formulation [l3| : 


(44) 


■mn 


V 

mm 

^mm 


r 

mm 

T 2 

^mm 


From Eqs. (43) and (44): 


^ J'/f • ip dv 
^mm ^ 

where ^ is the current density of source m and ^ is the free -space 
electric field. The self impedance of the longitudinal strip source 
(Eq. (20)), the transverse strip dipole f9] (Eq. (30) and Eq. (31)), 
and the rectangular surface dipole (Eqs. (39) and (40)), as a func- 
tion of size, are listed in Table I, II, and III, respecti vely. 

The mutual impedance between two sources is defined by 


(46) Z^n 


-1 
T i 
^mm *nn 


K 


En 


in dv. 


Tables IV, V and VI list the mutual impedance between two longitudinal 
strip sources, transverse strip sources [9j , and rectangular surface 
dipoles, respectively. 
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TABLE I 


Self Impedance of Longitudinal Strip Source 
shown in Fig. 9. 


h/x 

^11 


0.1 

582.38 

771.74 

0.2 

554.63 

488.50 

0.3 

511.90 

317.22 

0.4 

458.70 

199.36 

0.5 

400.65 

118.30 


CO 


♦ 



CD 


Fig. 9— Strip Source with £ = z . 
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TABLE II 


Self Impedance of Center-Fed Strip-Dipole 
shown in Fig. 10 


Sj = tj = h 



h/X = 0.05 

h/X = 0.10 

h/A = 0.15 

h/A = 0.20 

45^ 

0.11 -j 14.1 

0.47 -j 13.4 

1.20 -j 12.1 

2.53 -j 10.4 

90^ 

0.38 -j 20.9 

1.59 -j 19.3 

3.94 -j 17.1 

8.10 -j 14,0 

las'" 

0.64 -j 24.3 

2.68 -j 22.2 

6.49 -j 19.5 

12.95 -j 16.4 

I— » 

00 

o 

O 

0.75 -j 25.3 

3.11 -j 23.0 
1 

7.50 -j 20.3 

14.77 -j 17.4 


Fig. 10— Nonpl anar strip dipole with edges at Sj and t 
and terminals at 0. 











TABLE III 


Self Impedance of Center^Fed Planar Surface -Dipole 
shown in Fig. 11 


a/A 

Rii 


0.2 

11.48 

-69.76 

0,3 

23.98 

-35.26 

0.4 

38.80 

-15.36 

0.5 

52.98 

- 4.52 



Fig. U--Surface Dipole with ^ = z cos(iTz/a) 
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TABLE IV 


Mutual Impedance of Coplanar Strip Sources 
shown in Fig. 12 

hi = h2 = 0.1/X 


d/A 

^12 

Xl2 

0.0 

582.38 

771.74 

0.1 

526.74 

205.89 

0.2 

376 .25 

-140.97 

0.3 

171.77 

-283.09 

0.4 

-27.43 

-283.44 

0.5 

-172.11 

-190.34 



Fig. 12 — Coupled Coplanar Strip Sources with ^ ^ = 
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TABLE V 


Mutual Impedance of Center-Fed Planar Strip-Dipoles 
shown in Fig. 13 

Segment length: h/x =0.1 

Distance between midpoints: p/x = 0,3 


<(> 

3 = 0° 

3 = 45° 

3 = 90° 

6 = 135° 

Qo 

1.94 +j 0.81 

1.36 +j 0.39 

P 

o 

o 

+ 

o 

* 

o 

o 

-1.36 -j 0.39 

o 

o 

1.42 -j 0.58 

1.62 +j 1.07 

0.87 +j 1.80 

-0.38 +j 1.73 

60° 

0.39 -j 2.68 

0.90 -j 0.66 

0.87 +j 1.80 

0.34 +j 3.03 

O 

o 

CD 

-.13 -j 3.53 

-.08 -j 2.57 

o 

o 

o 

+ 

o 

o 

o 

0.08 +j 2.57 
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TABLE VI 


Mutual Impedance of Center-Fed Coplanar Surface-Dipoles 

shown in Fig, 14 

a/A = 0.5 b/A = 0.25 


0.75 

1.076 -j 7.913 

-2.925 -j 5.950 

-5.968 +j 1.761 

0,8871 +j 7.019 

0.50 

24.43 +j 5.997 

10.04 -j 9.321 

-9.568 -j 6.65 

-7,194 +j 8,971 

0.25 

53.10 +j 55.87 

29.54 -j 9.070 

-9.468 -j 19.33 

-15.74 +j 6.426 

0.0 

57.09 +j 13.23 

39.07 -j 22.49 

-8.659 -j 26.75 

-19.32 +j 4.263 


0.0 

0*25 

— 

0.50 

0.75 



! 

j^b 


Fig. 14— Coupled surface dipoles 
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CHAPTER V 

THE EXCITATION COLUMN 


The complex quantities in Eq. (9) form the excitation colunn 
in the matrix equation C„^ Ip = An,. Physically, A^ is the reaction 
between the impressed source and the test source m. These reactions 
are independent of the surface impedance or the dielectric coating. 
For arbitrary impressed sources which generate E,- in free space, 

Eq. (11) gives 


(47) 


Am = 




dv . 


The above expressions require numerical integration over test 
source m. 


A. PLANE WAVE ILLUMINATION (TWO-DIMENSIONAL TM CASE) 

If an infinitely long, electric line source, parallel with the z 
axis, is located at a great distance from the conducting cylinder, the 
incident field (B,- ,H-j ) may be regarded as a plane wave with 


(48) 


, jk(x cos + y sin ) 
Li = z Eq e 


where ^ is the angular coordinate of the source, and Eq is the in- 
cident electric field intensity at the origin. Fig. 15 illustrates an 
incident plane wave illuminating a longitudinal strip source with unit 
current density in the z direction. The integration in Eq. (47) is 
readily performed to yield 


(49) 



^ gJ’Pl 

jk cos (a- <j)j ) 


where = k(x.cos<j>^ + yjsin<j5^) and a is the angle between the positive 
X axis and the'^ vector directed from point 1 to point 2. 


B. PLANE WAVE ILLUMINATION (TWO-DIMENSIONAL TE CASE) 

Consider a magnetic line source, parallel with the z axis andlocated 
at a great distance away, illuminating a strip dipole with sinusoidal 
electric current density (Eqs. (30) and (31)) flow in the direction 
from 1 to 2 and from 2 to 3 as shown in Fig. 16. The incident electric 
field may be regarded as a plane wave with 


(50) 


i-i = - 




jk(x cos (h* + y sin ) 
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Fig. 15 — A plane wave illuminates a strip source 


In this case, Eq. ( 47 ) is readily evaluated to yield 

[e'^’^i-(cos kh.-j cos(ai~ <j). )sin 

(51) = -n Ho i 1 ! 

k sin khj sin(ai- ) 

[e’^'^-(cos kh2-j cos (a^“ 45 )sin kh^)e’^^^j 

+]-i *” ^ 

k sin kh2 sin(a2~4>^-) 


where hi and hz are the 'dipole segment lengths. The angle between the 
positive X axis and the vector directed to the terminals from point 1 
is denoted a^. Similarly a 2 is the angle of the vector directed to the 
terminals from point 3. 

C. PLANE WAVE ILLUMINATION (THREE-DIMENSIONAL CASE) 

If an impressed source is located at a great distance away from 
a surface dipole with current distribution z' cosUz'/a) as shown 
in Fig. 17, the incident field may be regarded as a plane wave with 

^ jk(x'sin 0 ; cos + y'sin el sin 4 ' + z'cos 0 I) 

E. = ^ e T 1 T i 1 . 
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Fig. 16 — A plane wave illuminates an electric strip 
dipole. 


From Eqs. (52) and (47), 


(53) 


. z')2iiab 


sin(X-,‘) cos(Y-,*/2) 

^(Yi2^ir2) 


where = 0.5 kb sinej sin<j>^ 
Y^* = ka cos q! . 
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Fig. 17 — A plane wave illuminates an electric 

surface dipole. 
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CHAPTER VI 
APPLICATIONS 


In the preceding chapters, attention has been directed to the 
general theory of the reaction integral formulation for scattering 
problems. Suitable bases and test sources have been defined and the 
mutual impedance between them has been evaluated. In the following, 
scattering from dielectric-coated cylinders with transverse magnetic 
incident wave will first be investigated. 

A. CYLINDERS WITH THIN DIELECTRIC COATING (TM CASE) 

Consider a dielectric-coated conducting cylinder illuminated by 
an electric line source. Let denote the surface -current 

density induced on the conducting surface, and denotes the 
polarization-current density induced in the dielectric layer. The 
first step in the reaction-Galerkin approach is to approximate the 
cylinder by a polygon cylinder with N segments. This is accomplished 
by fitting the cylinder with segments such that the perimeter of the 
polygon cylinder is equal to that of the original cylinder. 

Fig. 18a illustrates a dielectric-coated, conducting polygon 
cylinder illuminated by a parallel electric line source Ji z. Let 
Ip denote the current density ^ on segment n of the polygon cylinder. 
Each modal current Jn has a uniform current distribution as in Eq. 
(20). Now one represents ^ as the superposition of the N modal 
currents with weighting I^. This gives a piecewise uniform expansion 
for ^ with N unknown constants. The expansion for will be 
considered next. 

For a transverse magnetic source such as the electric line source 
J-j z shown in Fig. 18a, the electric field has only a z-component. 

For a conducting cylinder with thin dielectric coating. Maxwell's 
equations and the boundary conditions can be employed to obtain a 
suitable approximation for the electric field in the dielectric layer 
(Append! x)r 

-ki sin(kic)^ 

(54) £ = (TM case) 

jwe 


where ^ is the z-directed surface- current distribution on the con- 
ducting surface and e is the permittivity of the layer. Distance 
measured normally outward from the conducting surface is denoted by 
the coordinate c, and kj is the propagation constant in the dielectric 
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Fig. 18a- -Dielectric-coated, conducting polygon cylinder 
illuminated by a parallel electric line source. 


From Eq. (12) and (54), the polarization -current density can be 
expressed in terms of the surface -current density as follows. 

(55) *^1 sin(kic)^ (TM case) 

Since the surface -current distribution ^ has been expanded with 
rectangular-pulse bases one obtains a dependent expansion for the 
polarization -current density 

N 

(56) = I In 

n=l 
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where 


(57) 2i = s1n(kjc) J|| (TM case) 

The magnetic current PL vanishes if the cylinder is a perfect 
conductor. If the cylinder has finite conductivity and the impedance 
boundary condition (Eq. (6)) is employed, then the magnetic surface- 
current density can be expanded as 

( 58 ) M3 = ? I, 

n=l 

where 


(59) Mn ■ Zg ^ X n (TM case) 

and n is the unit normal vector directed into the source region. 

From the above discussions, the original problem illustrated in 
Fig. 18a can be replaced by its equivalent electromagnetic model, in 
that the dielectric-coated polygon cylinder is represented by an array 
of N hypothetical sources radiating in free space as shown in Fig. 18b. 
Each source is a modal current distribution S^, which is a collection 
of mode currents 2n and Mn. with weighting Ip. 



Fig- 18b^-Electromagnetic model of the problem shown 
in Fig. 18a. 
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Based on the zero-reaction concept discussed in Chapter II, if 
one places an electric test source in the interior region, as illus- 
trated in Fig. 18b, this interior test source has zero reaction with 
the other sources. To determine N current sanples, one makes N 
independent reaction tests. This procedure generates a system of N 
simultaneous linear equations. Numerical solution of this system yields 
a stationary result [17] for the samples of the current distribution. 

Let Cj^ denote the reaction between test-source m and mode 
current ^ in Fig. 18b. The reaction between the test-source m and 
the impressed source ;J-i plus the reaction between the test-source m 
and all the mode currents, N 

I ^ 

n=l ^ 

must vanish, leading again to Eq. (9) in which the reaction is 
given by Eq. (10) plus an additional term given by Eq. ( 14 ). 

B. CYLINDERS WITH THIN DIELECTRIC COATING (TE CASE) 

Consider a dielectric- coated, conducting polygon cylinder 
illuminated by a parallel magnetic line source I, as shown in 

Fig. 19a. The iiduced surface -current density flows in the di- 

rection transverse to z axis. Let I denote the current density at 
the corners of the polygon, and let one define N strip dipole -mode 
currents on the conducting surface. Mode 1 extends from point R to * “ 

point 2, Mode 2 extends from point 1 to point 3, and so on. Each 
mode ^ has a sinusoidal current distribution and unit terminal 
current density defined by Eqs. (30) and (31). Now one represents 
^ as the superposition of the N overlapping dipole-mode currents with 
weighting 1^, This gives a piecewise-sinusoidal expansion for ^ with 
N unknown constants. 

For the TE case, the electric field in the thin dielectric coat- 
ing is essentially normal to the conducting surface, and can be 
determined from the charge density distributed on the conducting 
surface. Via Maxwell's equations and the boundary conditions, a 
suitable approximation for the electric field in the dielectric region 
(Appendix) is 


(60) E = ~ cos(k ?) (z X j;) (TE CASE) 

~ ju)e ^ 

where k = k/e^.-! and b^, is the relative permi tti vi ty. Coordinate 
denotes^the distance measured normally outward from the conducting 
surface. ^ is the derivative of the surface-current density. From 
Eqs. (12) and (60), one obtains a suitable expansion for the polari- 
zation-current density inside the thin dielectric layer coated on a 
conducting cylinder as follows; 
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wi th 


(64) 4 = Zsfii ^ 


(TE case) 


Thus, the electromagnetic model, in the TE case, for the problem' 
shown in Fig. 19a is an array of N overlapping, si nusoi dal -dipole 
sources radiating in free space as illustrated in Fig. 19b, Each 
source is a modal current distribution which is comprised of 
the mode currents 5, and with weighting 1^. 



Fig. 19b--Electromagnetic model for the problem of 
Fig. 19a. 


Following the same argument used in the previous section, one 
performs N independent reaction tests by moving an electric test 
source to the conducting surface all illustrated in Fig, 19b, Again, 
the zero- re act ion concept leads to the matrix equation, Eq. (9), 
which can be solved nunerically for the stationary current samples. 
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In the reaction formulation, a polygon cylinder with N segments 
is employed to represent a cylinder with arbitrary contour. It is 
found that in order to obtain accurate results one must use at least 
five segments per wavelength, 

C, PLATES AND CORNER REFLECTORS WITH PERFECT CONDUCTIVITY 

Electromagnetic modeling of dielectric- coated, conducting bodies 
with finite extent can be accomplished through the procedure described 
in the previous sections. In this section, however, two specific 
cases, namely, scattering from perfectly-conducting plates and comer 
reflectors are discussed. The plates and comer reflectors are 
assumed to have an infinitesimal thickness. Therefore, the total 
surface-current density ^ (the vector sum of the current density on 
the front and back of the plates) is employed in the analysis. 

Consider the problems of plane wave scattering by rectangular 
plates and comer reflectors. The planar surfaces are divided into 
rectangular cells as illustrated in Fig. 20, TTie surface -current 
density is then expanded into two orthogonal sets of overlapping 
dipole -mode currents. Each dipole-mode current covers two cells 
and has a sinusoidal current distribution as defined by Eqs. (39) and 
(40). In Fig. 20 the arrows represent the mode current densities 
Thus, in the reaction calculation, the plates or the comer reflectors 
are represented by an array of overlapping mode currents radiating 
in free space and the reaction tests are enforced with a set of 
electric test sources. 

The zero-reaction concept leads again to the matrix equation, 

Eq. (9), which is then solved for the samples of the surface -current 
density. 
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Fig. 20— Electromagnetic model of plate and corner 
reflector. 
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CHAPTER VU 

FAR-FIELD RADIATION AND SCATTERING 


The field scattered by a conducting body is the sum of the free 
space fields generated by the electric surface -current distribution 
Jg and the magnetic surface-currents ^ induced on the conducting 
surface. If the body has a dielectric coating, the contribution from 
the polarization-current density should be included. To obtain 
the total field one adds the free space field of the impressed 
current source (2-,- ) . 

In the reaction calculation the continuous conducting surface is 
segmentized and the current distribution is represented by a set of 
modal currents. Thus, the field scattered by a conducting body is 
the sum of the free space fields generated by these modal currents. 
The far-field contribution due to three types of modal currents 
employed in this dissertation will be examined in the following 
sections. 

A. LONGITUDINAL, DIELECTRIC- COATED STRIP SOURCE 

For an electric line source of infinite length located on 
the z axis, the free space field is 

(65) E = -zkn Jz H^‘^(kp)/4 

(66) H = - ^ Hf’(kp)/4k. 


If the electric line source Is parallel with the z axis and passes 
through the point (x,y), its free space field at a distant point 
is 

■ — jkp 

^ kn dz^2j e jk(x cos <i> - y sin e) 

(67) E = -z e 

4/rTkp 

Fij. 21 illustrates an electric surface-current distribution 
J{t) = z J(t) on a dielectric-coated, conducting strip extending from 
rxi5yi) to (X 2 ,y 2 )* Distance along the strip is measured by the co- 
ordinate t, and distance perpendicular to the strip is measured by the 
coordinate c. For an arbitrary point on the strip. 


(68) X = Xi + t cos a 
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Fig, 21“-A dielectric-coated strip source. 
(^9) y = yj + t sin a . 


From Eq. (67), the free space field of this source ^(t) = z J^it) at 
a distant point (p, (f>) is ^ 


(70) 






dt 


where 


(71) = k(x^ cos (() + Yj sin<^) 

(72) c = cos (a- <j)) 


Using the approximation (Eq. (55)), the far-field contribution 
from the volume current density is 
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( 73 ) 



kr[/2j e 




4/7Tkp 


[ c " I 
V r / 


h ,*d 


ij 


0 -'0 


sin(kic) 


-jkcsin{a-4)) 
e dc dt . 


The integration along ^ can be readily performed to yield 


(74) 


^eq 


^TO - 


where 

(75) 






j[ki-k sin(a-(())]d -jlk^^+k sin(a-(i))]d 

e -1 e ■ -1 

^ + : ^ 


ki-k sin(a-(|)) 


k.+k sin (a- 45 ) 


If one introduces finite conductivity to the conducting strip and 
uses the impedance boundary condition (Eq. (6)), the free space field 
due to this magnetic current distribution is 


(76) 


eM = 

— n “ 


From Eqs » (70), (74), and (76), the distant scattered field 
generated by the dielectric-coated conducting strip (shown in Fig. 21) 
with 02 (t) = 1 is 


(77) - 

/Sirkp 


1 + 




z- sin(a-(|)) 


TM 


where 


(78) F 


TM 


_ e 


, e^>i 


(79) (|/2 = k(x 2 cos (fi + Y 2 sin(f)) . 

For broadside direction, i.e., a- 4 = ir/2 or 3it/2, 
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(80) = jkh . 


B. TRANSVERSE, DIELECTRIC-COATED STRIP MONOPOLE 

Consider again the dielectric-coated strip monopole shown in 
Fig. 21. This time the source is a transverse electric surface- 
current distribution 0(t) = t Jt(t). From Eq. (24), the magnetic 
field at a distant point (p,c|)) is 


(81) 


mJ _ A sin(a-<})) 
H — — z ——————— 

4Ak^ 


-jkp ji/j, ,h 
e e J^(t) 


jkct 
e dt 


From Eq. (62), the far-field contribution from the volume current 
density is 


/no^ cos(a-<}>) -jkp (c^^l 

(82) H = z e e ^ \~ — 

4Akp V . 


.h d 


a jkct -jk sin(a- (j>)c 

cos(k c) d!(t) e e ' dc dt . 

- 0 Q Z 


The integration along c is performed to yield 


zkv^ 


!i “ ^ cos(a-(})) e e 1 j!(t) e 

4/irko ^ n ^ 


4 /lilcp 


jkct 


dt 


where 

(84) 


TE 


f'er”!' 

I2j£y,, 


ej[k^-k s1n(a-4,)]d _j sin(c<- 


k^-k sin(a-(|)) 




k^+k sin (a- 


Consider a magnetic surface-current distribution M 7 (t) on the 
planar strip extending from (xivyi) to (x 2 ,y 2 ) as in Fig. 21. From 
duality relations in electric and magnetic systems and Eq. (70), the 
free space field of this source at a distant point (p, <{) is 


(85) 





4n A kp 


rh jkct 

Mz(t) e 

Iq 


dt . 
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If the magnetic surface -current density M,(t) in Eq. (85) arises 
through the finite conductivity of the cylinder, the impedance 
boundary condition in Eq. (6) yields 


(86) M^Ct) = s Zs Jt(t) 


where 


(87) s = (t X - 2 = ±1 


and the unit normal n is directed into the source region. 

From Eqs. (81), (83), (85), and (86), the distant scattered field 
from one segment of the dielectric-coated conducting cylinder (the 
segment in Fig. 21) is 

/ t c i^[(n sin(a-(j))+s Ze)Fi-n cos(a- (J>)gjr Foje'^*^^ 

( 88 ) = - — — — 

4n sin kh/irkp 

where 

/h jkct 

(89) Fj = k sin(kh) e J^(t) e dt 


and 


, , fb jkct 

(90) F 2 = k sin(kh) e ^ J|(t) e dt . 

If the transverse electric current on this segment has a 
sinusoidal distribution as follows 

sin(kh-kt) + K sin(kt) 

(91) J (t) = ^ , 

sin(kh) 

Eqs. (89) and (90) yield 
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(92) Fj = i— 

sin^(a-^) L 


-(cos kh + j c sin kh) 


I- 


and 


(93) = 


sin 2 (a- (j)) 

h 

sin^(a-{|)) 

I. 


e ^*^2 -(cos kh - j c sin kh) e '^’^2 






j c e 2 + (sin kh - j c cos kh) e ^ 


J>o 

j c e ^ - (sin kh + j c cos kh) e ^ 


sin2(a-(j)) 

In the end-fi re direction where (a-<t>) is zero or tt 


(94) Fi = 


■ J>1 . .. J>2 1 

e sin kh - kh e 

“J>2 J>1 

e sin kh - kh e 


j c Ij/2 
j c I 2/2 


and 


(95) F^ = 


kh e + sin kh e 


3*^1 . J’l '2 

kh e + sin kh e 


1 1/2 

1 2/2 


C. RECTANGULAR SURFACE DIPOLE 


^Consider an electric surface dipole with current density 
z' cos(7Tz7a) located on the y'z' plane as shown in Fig. 22. From 
reciprocity, the free space electric field generated by this source 
at a distant point (r, ^,.|)') (from Eq. (53)) is 


(96) 



‘'i 



a b sin 0 ' 
s 


sin(Xg) 


cos(Y^/2) g-jkr 

(YgS- n2) 


where 
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z 



0 

/ — ^ 

X 

Fig. 22 — A surface dipole radiates in free space. 


(97) Xg = 0.5 kb sin 0 ^ sin<^ 

(98) Yj = ka cos ^ . 

The e> and components of the scattered field with respect to the 
reference coordinate system 0 can be obtained easily via appropriate 
coordinate transformation. 

In plane-wave scattering problems, one is usually interested in 
the echo information defined as follows: 
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2 


(99) 




(TM case) , 

echo width = 

1 i m 27Tp ' 

p->«) 

1 eJ 1" 
|h"i^ 

(100) 

echo width = 

1 i m 27Tp - 

p-w 

IhM'' 

(TE case) , 


( 101 ) 


echo area = Tim 47rr^ 



(three-dimensional case)* 


In three-dimensional antenna problems, one is interested in the 
power gain: 


(102) gain 


W |E|2 
n| V|2 g 


where V is the terminal voltage and 6 is the conductance of the 
antenna. 
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CHAPTER VIII 
NUMERICAL RESULTS 


A. TM DIELECTRIC-COATED CYLINDERS 

In this section numerical results are presented for the back- 
scattering echo width of cylinders with thin dielectric coating. The 
incident wave is a transverse-magneti c plane wave. Fig. 23 presents 
the backscattering echo width of a circular cylinder with a thin 
dielectric coating. In the reaction calculation the cylinder is 
divided into N segments with N = 12 + 20 d/x and d is the diameter. 
Similar results for a square cylinder with 16 segments are shown in 
Fig. 24. For comparison, similar results are presented in Figs. 25 
and 26 for un coated cylinders with finite conductivity. Bi static 
echo width as a function of aspect angle is presented in Fig, 27 
through 29 for circular and square cylinders with a thin dielectric 
coating. In all cases, the dielectric layer has a relative dielectric 
constant of 10. 

B. TE DIELECTRIC-COATED CYLINDERS 

Numerical results ai^ presented for cylinders with a thin 
dielectric coating illuminated by a transverse -electric incident 
plane wave. Fig. 30 presents the backscattering echo width of a 
circular cylinder with a thin dielectric coating. As before, the 
cylinder is divided into N segments with N = 12 + 20 d/X and d is the 
diameter. Similar results for a square cylinder with 16 segments are 
presented in Fig. 31. In Figs. 30 and 31, backscattering echo width 
for uncoated cylinders are also included for comparison. Bi static 
echo width as a function of observation angle for coated cylinders 

are shown in Figs. 32, 33 and 34. The relative dielectric constant is 
1.5 in each case. 

Examining the results, one can note that the reaction calculation 
gives accurate data for uncoated conducting cylinders. Satisfactory 
results are also obtained for cylinders with a thin dielectric coating. 

C. PERFECTLY-CONDUCTING PLATES AND CORNER REFLECTORS 

Fig. 35 presents the backscattering echo area of a square plate 
with perfect conductivity for the broadside aspect. In the reaction 
calculation, the plate is divided into cells, and overlapping current 
modes were employed as illustrated in Fig. 36. In this case the 
transverse current was neglected and 45 modes were used for the current 
distribution. Useful results can be obtained with as few as one mode 
per square wavelength of surface area. For comparison. Fig. 35 also 
shows the experimental measurements of Kouyoumjian [18]. 

The magnitude and phase of* the induced current density on a 
perfectly-conducting rectangular plate are illustrated in Figs. 37 and 


45 



ECHO WIDTH / WAVELENGTH 


1.4 



Fig. 23— Backscattering echo width of dielect r1 c-coated 
ci rcul ar cylinder , 
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Fig. 26- -B roads Ide backscatteri og echo width of square cylinder. 
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Fig. 2 . 1 — Bi static echo width of dielectric-coated circular 
cylinder. 
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Fig. 29 -Bistatic echo width of dielectric-coated square cylinder for oblique incidence. 
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Fig. 30--Backscattering echo width of dielectri c-coated 
ci rcular cyl inder. 
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Fig. 31— Broadside backscatteri ng echo width of dielectric- 
coated square cylinder. 
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Fig. 32- -Bi static echo width of dielectric-coated circular 
cylinder. 
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Fig, 33--8ist.atic echo width of dielectric-coated square 
cylinder for broadside incidence. 
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Fig. 34--B1stati c echo width of dielectric-coated square cylinder for oblique incidence. 







38. Figs. 39 through 42 show the normalized backscatter cross 
section of a rectangular plate. Figs. 43 and 44 show the normalized 
backscatter cross section of a comer reflector. The title of each 
figure gives the echo area at the broadside aspect in terms of dB = 
10 log(a/x2). 

Figs. 46 through 49 show the E-plane gain of the comer- 
reflector antenna illustrated in Fig. 45. Figs. 50 through 53 show 
the H -plane gain of the same antenna. For comparison. Figs. 46 
through 53 include experimental measurements obtained by Melvin 
Gilreath at NASA Langley Research Center. In the experimental 
measurements the receiving antenna was linearly polarized in the 
theta direction. Similarly, the calculated gain is based on E . 

The dipole length is x/2 and the radius is 0.005X. 

In the reaction calculation, only vertical modes were employed 
to approximate the current distribution. The number of modes used 
to obtain the results given in Figs. 37 through 53 are listed below. 
In each case, the matrix size is equal to the number of modes. 

Fi gs . Number of Modes 


37,38 45 

39,40 55 

41,42 75 

43,44 30 

46-53 61 



y/x 


Fig. 37 — Magnitude of surface-current density induced on 
a perfectly-conducting rectangular plate for a 
plane wave incident at broads ide. 
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Fig. 38-“Phase of surface -current density induced on 
a rectangular plate for a plane wave incident 
at broadside. 



yz plane of a rectangular plate, a , J e, i) ^ 
15.25 dB at (90^,90°). 
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Fig. 40-'-Nonnalizecl backscatter cross-section 1n the 
yz plane of a rectangular plate, ooi e, d = 
15.23 dB at (90^,90''). 



Fig. 41— Normalized backscatter cross-section in the 
xy plane of a rectangular plate, a.ie, a) = 
15.25 dB at (90°,90“). 
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Fig. 42— Normalized backscatter cross-section in the 
xy plane of a rectangular plate. aonCe, 6) 
15.23 dB at (90°, 900). 
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Fig. 43— Normalized backscatter cross-section in the 
xy plane of a comer reflector, 0 , d,) = 

-0.38 dB at (90'', 90^). 



Fig. 44-’-Normalized backs catter cross -section in the 
yz plane of a comer reflector. aoo( 0 ,(b) = 
-.38 dB at (90^90°). 0 8"^ ‘f' 
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Fig. 45 ““Comer- re fleeter antenna. 
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Fig. 46— Relative gain in the E-plane of a comer- 
reflector antenna. G{ e, (jj) = 4.31 dB at 
(90°, 90°). 
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• • •REACTION S=0.75X a^lSO” 



Fig. 48 — Relative gain in the E-plane of a comer- 
reflector antenna. G( e, <b) = 7.48 dB at 
{ 90 ^, 90 ^), 
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Fig, 49— Relative gain in the E-plane of a corner- 
reflector antenna. G( e, 4) = -1.06 dB at 
(90^,90^). 
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Fig. 50“-Re1ative gain in the H-plane of a corner- 
reflector antenna. G( e, 4 = 4.31 dB at 
(90^y90^). 
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•••REACTION S =0.75X, a =180° 



Fig. 52— Relative gain in the H-plane of a comer- 
reflector antenna. G( e, <j)) - 7.48 dB at 
(90°, 90°). 
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Fig. 53—Relative gain in the H-plane of a comer- 
reflector antenna. G( 8, (|)) = -1.06 dB at 
(90°, 90^). 


72 


CHAPTER I X 

SUWARY AND DISCUSSIONS 


The reaction concept and Galerkin's method are employed to 
develop an integral equation formulation for radiation and scattering 
from perfectly-conducting plates, comer reflectors and dielectric- 
coated conducting cylinders. 

For the two-dimensional problems, the contour of the cylinder 
is divided into segments and the surface -current density on the 
conducting surface is expanded with sinusoidal bases for the TE 
polarization and rectangular-pulse bases for the TM polarization. 
Reaction tests are enforced with electric test sources. This reaction 
Galerkin technique yields accurate results for scattering by cylinders 
with as few as five segments per wavelength. Furthermore, this 
technique provides a symmetric impedance matrix in the matrix equation 
The point-matching procedure, on the other hand, generates an unsym- 
metric impedance matrix and requires on the order of ten segments per 
wavelength to yield accurate results. 

For a coated cylinder, the dielectric layer is modeled with the 
equivalent polarization current radiating in free space. Maxwell’s 
equations and the boundary conditions are employed to express the 
polarization-current distribution in terms of the surface-current 
density on the conducting surface. The impedance matrix has the same 
size for an uncoated cylinder and a cylinder with a thin dielectric 
coating. It is found that the polarization-current model is more 
accurate than the popular surface -impedance model. 

For the three-dimensional problems, a new model that involves 
dividing the conducting surface into cells, expanding the current 
distribution with subsectional bases and enforcing reaction tests with 
electric surface dipoles is developed and applied to the problems of 
radiation and scattering from perfectly-conducting rectangular plates 
and corner reflectors up to six square wavelengths in size. For arbi- 
trary aspect and polarization, this sinusoidal-Galerkin technique 
yields good results for scattering by a one -wave length square plate 
with an 18 x 18 matrix. For this same problem, the wire-grid model 
(and the surface-cell model with pulse bases and collocation) requires 
a matrix size of 60 x 60. Thus, the new surface-current model offers 
a substantial improvement in computer storage requirements. This 
model can also be applied to nonplanar surfaces provided that the 
general surface is approximated by a set of planar cells. In this 
case, the cells may have to be reduced in size to obtain a satis- 
factory fit. This will increase the number of cells enployed for 
modeling, thereby limiting the application of this model for non- 
planar surfaces. However, extensions can be made such that the indi- 
vidual cell is doubly-curved in nature, and integration techniques can 
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be developed for evaluating the reaction between these curved cells. 
With this extension, the problem of radiation and scattering from 
curved surfaces can be analyzed more efficiently, and bodies with 
surface ai^as up to ten square wavelengths may be analyzed. For 
problems involving more complicated geometry, such as antennas 
mounted on aircraft, a combination of planar and curved cells can be 
employed to model the aircraft surfaces. 
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APPENDIX A 


ELECTRIC FIELD INDUCED IN THE THIN DIELECTRIC LAYER COATED ON A 
PERFECTLY-CONDUCTING POLYGON CYLINDER ILLUMINATED BY AN INCIDENT 
PLANE WAVE 

Consider a dielectric-coated, perfectly-conducting polygon 
cylinder illuminated by an incident plane wave as shown below. The 
dielectric layer is a source-free region and has a thickness of d 
and a dielectric constant e. The segment length of the cylinder is 
denoted by i. Define a coordinate system such that Z x t = z, 



Fig. 54 — Dielectric-coated, perfectly-conducting polygon 
cylinder illuminated by an incident plane wave. 


where t is a unit vector tangent to the conducting surface and Z is a 
unit vector normal to the conducting surface. The incident plane wave 
is either transverse-electric or transverse-magnetic with respect to 
t-ax1s. 

For a perfectly-conducting cylinder with a thin dielectric 
coating (d/t « 1), the magnetic field inside the dielectric layer 
can be expressed as follows: 
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(103) H=HsCOs(3c) 


where ^ is the magnetic field induced on the conducting surface S and 
is related to the surface -current density ^ by 

(104) X ^ 


Coordinate c measures the distance normally outward from the conducting 
surface and 6 is the transverse propagation constant in the dielectric 
region which can be determined from the wave equation. From Eqs. (103) 
and (104) and Maxwell's equations, the electric field inside the 
dielectric layer can be given as 

(105) E = X c) cos(bc)]. 


For the transverse-magnetic incidence case, the surface current 
density has only a z-component ^ = z Js(t) and e = = k/e^. From 

Eq. (105), it can be shown that 

(106) i sin(kic) i (TM case) 

For transverse-electric case, the surf ace -cur rent density has 
only a t-component ^ = t Js(t) and 6 = k = k/e^ V . For this case, 
Eq. (105) yields ? 

(107) £= x^j') cos(k^?) + sin(k^ 5 ) . 

If « Is Eq. (107) reduces to 

(108) E = cos(k^c) (z X ^') (TE case) 

where is the derivative of the surface-current density. 
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APPENDIX B 

COMPUTER PROGRAMS FOR RADIATION AND SCATTERING 
FROM TM DIELECTRIC-COATED CYLINDERS 


INCLUDE. CROUTB,2902W 
COMPLEX C(40,40) *ZS|CJ(40) 

COMPLEX FJ,PJJ,CST,EM,EJE 

DIMENSION XI^O) ,Y(40) tD(AO) , I A (40) ,IB(40) 

C0MM0N/C0A/CCC,FR2, TSK2L 

DATA IDMt INT/40t 10/ 

data PI , TP/3. 14159, 6 *28318/ 

WRITE! 1,998 I 

998 F0RMAT(5x,*L0P=?,2S=?/* ) 

RPAD(0,-) LDP,ZS 

10 Continue 

WRITE! 1,999) 

999 FORMAT! 5X, 'NCASE=? 1 = RFCTANGULAR,2=C IRCULAR , 3=RTR IP/ • ) 
READ(0,-) NCASE 

WRITE! 1,997) 

READ!0,-) PH1,DPH,BSC 
GO TO !400, 500, 300) , NCASE 
300 CDNTiNUE 
PHI=90. 

READ(0,-) WK,NM,TSL »ER2 
TSK2L=TSL«TP«SQRT(FR2> 

CCC=(EK2-1. )/ER2Y(l .-COS! TSK2L) ) 

WRITE !6 ,4) WK 
WK.= WK*2. 

NP=NM+1 

DX=WK/Nm 

DO 1 1 = 1, NP 

Ym = .o 

X!I ) =DX’M 1-1 ) 

1 CONTINUE 

DO ? J=1,NM 
IA!J)=J 
IB(J)=J+1 
0!J)=DX 

C J! J)=CMPLX!2 .♦PV/DX, .0) 

2 CONTINUE 
GO TO 600 

500 continue 

PHI=.0 

PHR=PHI^PI/180. 

CPH=COS!PHR) 

SPH=S1N!PHR) 

CST = TP/(-30.=«^! .707,-*707) ) 

READfO,-) DL,TSL,ER2,XKS 
WRltF<6,4) DL,TSL,FR2 
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cc 


TSK2L=TSL*TP^S0RT( ER2) 
CCC:t(ER2-“l*) /PR2^(1 .“C0S{TSK2L> ) 


CC 

NM=I2-^20.♦0L 

1F(MK.LT*16) NM=16 

WRIT£a»-)MR 

NfP = NM 

PHO=TP/NK 

Pl^DL/NM) /SIN( .5*PH0) 

DO 11 I^lrNM 

X( I )=TP*BL^CDS(PHO*(I-l) ) 

Y{I)=TP + BL^SI\’(PHO*( I-U ) 

IA(I)=1 

1B( 

IF(I.FQ.NM) IB{U=1 

D{ 1 ):rTP*FL^2.’>^S •5<'PHO) 

11 CONTIWUe 
GO TO 80 
AGO CONTI?^UC 

^^7 FOPMAK 5X, 'PHI ,nPH,BSC = ?,BS,C>0. ♦♦YES/') 

IF(L0P.EQ,2} READ(0,~) XKS,YKS 
I F( LOP <,£0.2 ) WR ITE (6fA)XKS ,YKS 
PHP=PHI *pi/ieo. 

CPH=COStPHR) 

SPH-SIN(PHr) 

CST=TP/(«30.*( .707,-.707) ) 

REA0(0^--) AXv PYvTSLpER2 
WRITF(6»A) AX»RY,T$Lf ER2 
T$K2L=TSL«Tp*$QRT(ER2) 

CCC=(Fr2-1. )/FP2’^(1--C0S(TSK2L) ) 

CALL REC ( AXfBYt IDM,^^MtNP , X ,Y , I A , IB ) 

DO 260 I=ltNM 
xn)=TP«x(i) 

Y(I)=TP*Y(I) 

260 CONTINUE 

DO 45 J=l,NM 
K=IA(J) 

L=1B(J) 

D(J)=SQRT( (X(U-X(K) )**2+(Y(U-Y(K))*«2) 

45 CONTINUF 
80 CONTINUF 

DO 2? 1=1, NM 
KA=lA(n 
KB=IBU ) 

IF<LDP.F0.2) GO TO 33 

CALL CFF(X(KA),Y(KA) ,X(KB),Y(KB) tOn ) f CPH,SPH,2S,EJ»FMfFJF) 
C J( n = FJ*CST/0( I) /D( I) 

CCC 

GO TO 22 
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33 CnNTINUE 

C &L L c FL S ( X ( K A ) t Y ( K A ) » X ( K B ) , Y ( K B ) , XK S f YK S t D ( n 1 1 N' T, F J ) 

C J( I ) = F J^TP#TP/D( n/0( I ) 

c j( n=cju )/(-6o.+pi) 
ccc 

WP1TF(6,4) Cjd) 

22 CONTINUE 

600 continue 

lSYf^=0 

C 

R$=CABS(2E) 

1F(RS,GT.0.) 1SYM=1 
C 

cc 

IF(TSL.C1,0*) ISYN?=1 

CC 

CALL CDANT(CtO,XrY,2S, IA,IB, ISYM, rOM, INT»NM,N?I 
1F( ISYM.FQ.IO) GO TO 1000 
DO 3 1=1,1 
DO ? J=ltNN 
WRITF(6,4) C ( I, J) 

4 FORMAT! 3X,BF10.4/I 

3 CONTINUE 

CALL CRaUTlC,CJ,NM,IDMtISYM,l,l) 

PM=0*0 

I FCBSC*GT*0*) PH=PHT 
30 CONTINUE 

PHR=PH4PI/160. 

CPH=C05IPHR) 

SPHsSIN(PHR) 

FJJ=(.0,»0) 

on 250 K=1,NM 
KA=IA(K) 

KP=IB(KI 

CALL CF F ( x(KA),Y(KA I, X( KB ) ,Y( KB) fD(K),CPH,SPH, 
?Z$tFJ,FM,FJF) 

F; JJ = FJJ+ ( EJ+F JE )«C J(K ) + FM^CJ (K) 

250 CONTINUF 

IF( LOP. FO.? ) FJJ = FJJ-*-CEXP (CMPLX ( .0 , XKS *C PH+YK S*S PH+ P I /4. ) ) 
EA3=CAe$( FJJ) 

FWL = 2. + PI’«'FAB*FAS 
EWS = tWL=*PI*2. 

WRITE! 1,-) PHtEWL 
WRITE (6, 4) PH,EWL,FAB 
ph=ph+dph 
PHFND=360. 

I F( BSC .GT.O. ) PHEND=PHI 
I F{PH*LE*PHENM go to 30 
1000 CONTINUE 

RFAD!0,~) IC 
IFnC.EQ.O) GO TO 10 
CALL r-X IT 
END 
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SUB^^OUTINF CDANT<C.D, Xf YtZS, I At 16, ISYM, IDM, INt,NM,NP) 

COMPLFX 2S,PU,C(IDM,1DI«») 

DIMENSION X( ir>M ) ,Y( IDM) ,0( IDH) , 1A( IDM) t IB nOM) 

00 20 1=1 ,NM 
DO 20 J=1,NM 
20 C(I,J)=^(.O,.0) 

0MAX = ,0 
no 23 J= 1 ,NM 
DK=0(J) 

25 IF(nK,GT*DMAX) OMAX=DK 
WPIT»^( 1 ,-)DMAX 
I F( DMAX *LT.3* ) GO TO 30 
ISY.w=10 
RETURN 
30 CONTINUF 

no 200 K=l,NM 
KA=IA(K) 

KB=IB(K) 

DK=n(K) 

LL=1 

IF( ISYM.FO.O) LL=K 
DO 200 L=LL,NH 
LA=IA(U 
LB=IB(L) 

DL=D(L) 

1F(K. FQ.U GO TO 120 
I ND=(LA~KA)*( LB-KA)*(LA~KB)*(LB~KB) 

1F( INO*EO.O) GO TO 80 

CALL ZMMCIXCKA) , Y ( K A ) ,X ( KB ) , Y ( K9 ) , X { LA ) , Y ILA ) ,X ( LB ) , Y (L B ) , 
^ZStOKtOLtlNTtPll) 

GO TO 168 
80 CONTINUE 
JM=KB 
JC=KA 

1ND=:<KB-LA) ( KB-“LB ) 

IF< IND.NF*0) GO TO 82 
JC=KR 

JM=KA 
82 JP=LA 

IF(LB.FO-JC) GO TO 83 
JP=LB 

8 3 Call ZMMBIXIJM) tYUMl ,X( jo ,Y( JCJ,X(JPI ,Y(JP),ZS,DK,DL, INT,P11) 

G 0 To 1 68 

120 CALL ZMMAIDK tZStPllJ 
168 C(K,L) = PU 
200 CONTINUE 
RETURN 
e ND 
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SUB POUT IMF 2MMA(DK,ZSfPl 1 ) 

COMPLEX 2S*HO,HltVH,GI*Y,PU 
CGMMON«/COA/crX,FR2TTSK2L 
DATA Pl/3. 14159/ 

TP=2^«PI 

FTA=120.’«^P1 

CALL HAMKIDK 1 ) 

GUVHIOK) 

D=0K/TP 

Pll=n/{ ./(2.^TP’5'TP)*( DKXcH 1-( .Gt2. ) /PI ) 

P11=P1 1’>^TP*FTA/D/D 

Y^Pll’^2./FTA 

cc 

Pll = Pll + Pil’>CCC 
cc 
c 

p11 = PU4ZS*PI/DK 

c 

PFTURM 

CMD 

SUBROUTIMF ZHPP ( X 1 , Y1 , X2 , Y2 t X3 » Y3 , ZS » 

?DK1,DKZ,1MT,Q11) 

COMPLEX VH,ZS,OU*Yll ♦ t CP , Gl , G2 , G12* 

2H 10,Mn TH20tH21,H120tH121»HPG,HPl,HW0,HMl 
COMPLEX DOllP ♦D0nMtPKHlP,p.KHlM,D011 tOCNT 
C0MM0N/C0A/CCC,EP2,TSK2L 
DATA CCP/( .0, .63662)/ 

DATA PI/3.14159/ 

PS=CABS (2S) 

CBET^I X2-XD/DK1 
SBET=(Y?-yi)/OKl 
XB=(X?-X1 )yC BE T4(Y3-Y1 ) ASSET 
YB=-(X3~X1 ) aSBET+(Y 3-Y1)*CBET 
C AL=(Xe -DK.l) / 0K.2 
S AL=ABS(Y6/0K2) 

AL=ATaN?< SAL^CAU 
CNT=15. A4.yPl AP I/DK2/DK1 
DK12=DK1+DK2 
lElCAL.LT.O.) GO TO 20 
IFC SAL .GT. 0.04) GO TO 20 
CALL HANK (OKI ,H10,H11 ,1 ) 

CALL HAMK(0K2,H?0»H2l,l ) 

DK12=DK1+DK2 

CALL HAMK(OK12fHl20tH121,l ) 

&1=VH(DK1) 

G2=VH(DK2) 

G12*VH( DKi2) 

Q 11=CNTa (~CCP**'DK1A (h 1 1-Gl ) +DK2A ( H2 1-G? ) ~DK1 2* ( HI 21-Gl 2) ) 

CC 

011=01 1+011 ACCC 

cc 
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c 

c 


c 


RETURN? 

20 CONTINUE 

iNp = ?»( iNr/2 ) 

I P= IMP+ 1 

JP-IP 

FIT=1NP 

ALT=AL/2. 

CALT=COS< ALT) 

SALT=SIN(ALT) 

RCP=nKl 2>C ALT 
RSP=- (DK2-DK1 ) ’►S ALT 
PHG=AT AN? (R SP,RCP) 

PHW=“ALT 

php-phc 

DPHM=(PHC^ALT>/PIT 
DPHc^ = ( ALT-PHO/FIT 
YU=(*0, »0) 

ncn i=r( . Ot.o) 

DO 200 1=1, IP 
D^SGI+3. 

1F( I .EQ*l.nR.I *F0*1P ) 0=U 
SAP = SIN ( ALT + PHP > 

5 AM=$1N< ALT-PHM ) 

APGP=DK24cSAL/SAP 

4PGM=0K1’5^SAL/SAM 

CALL HANK ( ARGP,HP0,HP1 , U 

call hank (ARGM,HM0,HM1, 1 ) 

1F( RS«LF.0« ) GO TO 300 

0APGP=APGP/F1 T 

DARGM=ARGM/FIT 

R KP = • 0 

RKM=*0 

SGJ=-1 . 

D0UP=( .0,.0) 

0011M=( .0, ,0) 

RKH1P=CCP 
RKH1M=CCP 
DO ICC J=1,JP 
C=SGJ+3. 

IF( J,F0*1*0R, J^FQ.JP) C=l. 
IF(J.E0.1) GO TO 9A 
CALL HANK, (RKP ,H10,H11,X ) 
CALL HANK(RKM,H20,H21il ) 
RKH1P=RKP*H11 
RKH1NI=RKM^H21 
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CONTINUE 

001 1P=DQ11P+RKH1P*C 
■ DailM=DQllM+RKHlM^C 
RKP-RK P-1'0 ARGP 
RKM=RKM+DARGN 
SGJ=:-SGJ 
100 CONTINUE 

DOl 1-DC 1 1 +S AM (001 1P*DARGP*DPHP + D011M’<'DARGM<‘DPHM ) 

300 CONTINUE 
C 

Y U=Yll + { ARGH*HMl»DPHM-»-ARGP«HPl^DPHP )*0 
PHM=PHM+DPHM 
PHP^PHP-i'DPHP 
SGI=:-S&I 
200 CONTINUE 

01l=CNT^{Yll/3./SAL-CCP’«‘AL/SAL) 

cc . 

QU=01 1+Qll^CCC 
CC 
C 

OCNT=ZS<fCMPLX( .Ot-2«*Pl)/ (4.^^DKl*r)K2*SAL) 

0Oll=DOn*DCNT/P. 

Qll“011+D011 

C 

RETURN 

ENO 

SUBROUTINE ZMMC(X1 ,Y1 ,X 2t Y 2, X3 t Y3 t X4, Y4, Z S t 
2DK1 fDK2 t INT,P11 ) 

complex Pll tG3,G13tC23#G123,VH »HAOf HA 1 ♦HBOtHS 1 ,HCO, HC If HD Of HD 1 , 

COMPLEX HXfHYfDPll fDCNT 

COMPLEX HO, HI 

COMMON/CDA/CCC,ER?,TSK2L 

DATA P1/3.1A159/ 

RS=:CABS(7S) 

CB«=T = (X2-X1)/0k1 

SBET=( Y2-Y1 ) /OKI 

XA=( X3-X.1 )*CBET + ( Y3“Y1) ’f^SBET 

XB=(X4-X1 )^CBFT-KY4~Y1)’^SBET 

YA=-( X3-xn=J‘SBET+(Y3-Yl )*CBFT 

YB=-(XA~X1)*SBET+(Y4-Y1)*CBET 

C AL=(XB-XA) /DK2 

SAL= (Y6-YA I/DK2 

CNT-15 •*4,+PI*PT/DK1/DK2 

ASAL = ABS( SALI 

I E( ASAL.GT. .04) GO TO 20 

IE(YB.NE..O) GO TO 20 

DK3=ABS(XA-DK1) 

c 


85 



c 


IF( XA.LT.0*0) 0K3=A8S(XB) 


cc 

cc 


DK13=DKl+DK3 

DK23=DK?+DK3 

DK123=0K1+DK2+DK3 

CALL HANK(nK3 »HAO,HAl,l) 

call HAMK(0K13 

CALI HANK(DK2? »HC0,HC1,1) 

CALL HANK(DK123 ,HDC,HD1 , n 
G3=VH(HK3 ) 

G13=VH(DK13) 

G?3=VH(nK23) 

G123=VH (DK123 ) 

Pll=CK’T«(DK3«iG3--HAl)~DK13’«‘(G13«HBl)-0K23^{G23~HCl)4'DK123* 

?(G123-HDD) 

P11=P11+P11*CCC 

PFTUPfVi 
20 CONTINUE 
RMIN=10000, 

X = XA 
Y =Y A 

DX=DK2+CAL/A. 

DY-DK2+SAL/4. 

DO 40 J-lf5 
Y<;=*o 

R=AeS(Y) ./ 

IF(R.GT .l.E-15 ) YS=Y*Y 

xs=o.o 

XARtrABS (X-DKl ) 

IF(XAB*GT*1*F-15) XS=XAB^XAB 
1F(X.LT.0,) R=SORT(X*X+VS ) 

IF(X.GT.DKI) R=SQRT{XS+YS) 

IF( R.LT.RMIN) RMINf-R 
X=X+DX 
Y=Y+DY 
40 COMTINUE 

FNT:=1* (4«INT) /lO 
ISS-FA’T^DKl/RMIM 
ISS=2+( JSS/2) 

1 F( ISS«LT.2) ISS=2 

IF(ISS.GT,20) 1S?; = 20 

FSS = ISS 

ISQ=ISS+1 

DS=DK1/FSS 

I TT = FMT’«^DK2/RI^IM 

1TT“2^( ITT/2 ) 

IFdTT.LT. 2) ITT^2 
IF< ITT.GT.20) lTT=r20 
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Frr= ITT 

1T0=^ITT + 1 
DT-DK 2 /FTT 
. 0 X=DT^CAL 
DY=DT*SAL 
X=XA 
Y=YA 
SGJ = “1 . 

PX 1 =(. 0 ,.D) 

C 

DPU = (. 0 ,* 0 ) 

C 

OD 200 J=^ 1 ,IT 0 
n=SGJ*»' 3 . 

I F( J. 60 * 1 «OR.J«EC.ITQ) 0 = 1 . 

XP =.0 

YS =*0 

YA 8 =ABS(Y) 

IF( YAB.GT. l.E- 15 ) YS=YA 6 *YA.B 
FZ=(. 0 ,* 0 ) 

C 

HX=< * 0 ,. 0 ) 

HY=(. 0 ,* 0 > 

C 

SG 1 =- 1 . 

DD 100 1 = 1 , ISQ 
C=SGl+ 3 » 

1 F( 1 .F 0*1 •OR. UEQ-ISQ) C = 1 , 
0 FLX=A 6 S(X“XP) 

nxs=*o 

IF(DFLX.GT* 1 ,F~ 15 ) DXS^DFLX’J'DFLX 
RK=SQRT(DXS+YS) 

C 

SPH=Y/RK 
C PH= ( X-XP ) /RK 
C 

CALL HANK(RK,HOtHl, 0 ) 

. EZ=EZ4H0+C 
XP=XP4DS 

SG 1 =~SGI 

C 

1 F(RS.LF. 0 ,) GO TO 100 
HX=HX+H1’J^(-EPH) *C 
HY=HY+H 1 «CPH*C 
C 

100 C0M71^^lfF 

EI=FZ^0S/3* 

Pll=PXl4BZ*D 

C 

HX=HX 4 DS/ 3 * 

HY=HY>»^DS/ 3 , 

DP 11 =DP 1 1 +(HX^CAL+HY^SAL )*D 
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SGJ^-SGJ 
X=X+DX 
Y=Y>DY 
20n CPi^JTlNUE 

Pil=Pll*nT/3. 

P11 = P11’^CMT 

c c 

P ll^Pl 1+Pll*ccc 
cc 
c 

DCtSJT-IS^CMPLX ( *0i-2»^P\ ) / ( ^.^OKl’^DK^ ) 

DPI 1 =0P1 1 ’>DC\*T*DT/3. 

P11=P1 1+OPll 
C 

RFTURM 

r-NO 

SUBROUTINE HfiNK{X»H,Hlt ID) 

COMPLEX H,H1 

DATA TSP/ .63661977/ 

1F( X.GT.3.)GD TO 100 
XLN=TSP^AL0G(X/2.) 

P = , 0 
B1=.0 

Y = . 0 
Y1=.0 
Xl=X/3. 

X?“Xl’i^Xl 

1 F( Xl.LT. .DGO TO 60 

X4=X2»X2 

X6=X2’’'X4 

lFtXl.LT..3)CrO TO 55 

XB=X2*X6 

X10=X2^X8 

X12 = X2’>X1G 

P=.21F—3^Xl?“.39444E-2^X10-^.444479E-l^XB 
Y=-*24BA6E-3*X12 + .427916E-2^X10-.4261214E-1=^XB 
B1=.1109F“4YX12-.31T61E-’3*X10+.443310F*‘2«X8 

Y 1= ,27B73E— 2’!'X12~.400976e-l’5'X104-#3123‘>5l=<‘X8 

55 B -B*“ • 3 163B66^X6 + 1 • 2 6562 08^X 4 
Y=Y>f .25 3001 17 YX6-. 743 503 84*X4 
B1=B1-# 3<554289E~1«'X6+.21093573*X4 

Y 3 = Y1~1 .3 l64827*X6-^2 * 16B2709 + X4 

60 B=b-2, 2499997«X2+1 . 

Y=Y + *605 59366’^X2+ .3 6746691 + XLN*B 

B1-X*( B1-.5624«>985*X2>.5 ) 

Y 1 = (Y1+ . 221 2091*^X2-. 6366 198 ) /X-^XLN^Bl 
GD TO 200 

100 SW=SQRT(X) 
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X1-3./X 

X3^X1^X2 

X42:X1^X3 

X5^X1*X4 

X^'=X1*X"' 

F:^^7 97eB456-‘*77F-6*Xl~.5 5274e~2=«'X2~,951?F-4*X3+. ISVeBTF-a’^'X^. 
2~*72805E-3^X5+, 14476F~3^X6 

T = X-^78 5398 16-*4l6f>397E-l«XI-*39 54E-4*X2^*262573F-2*X3 
2-.54125E~3*X4-.79333F-3*X5+. 13558E-3*X6 
3=F«CDS(T)/SW 
Y=F*SIN(T)/SW 

F -,7<^78 84b6+. 136E”-5^X1+. I6?>96f 7E~X*X2+ . 17103F-3*X3-*2 4951 1E-2*X^ 
2+ *113<S5 3E-2«X5-*20033F“3*X6 

T =X~2r3 56l9 43+.12499612«^Xl**’.565E-4’»^X2-.637879E-2*X3^* 7434eE-3^X^ 
2+.79824E-3<fX3-.29166E-3^X^ 

61 = F’^C0S(n/SW 
Y1=F^S1N( T) /SW 
200 H=CMPLX (B»-Y) 

Hl=CMPLX(Blt-'Yl) 

RETURN 

END 

COMPLEX FUNCTION VH(X) 

DIMENSION A(ft) tB(8) 

COMPLEX GI51 ) ,H0,H1 

data g/- 

?( .0, .0) , ( .19933,“.34E70) » < .39470»-.3093?) , 

? { .38 224,^.59927 )t ( .75834, -,63787 ), ( .91 973 , 63707 I , 

2( 1,063 5 6, -.60490) ,I 1.18750,^.54783) , ( 1 .2896 2 , -.47156 ) , 

2( 1.36940, -,381 36) , ( 1 .425 77, -.282 19) ,( 1 .459 13, -.17871 ) , 

2i 1.47029, -.07527) ,( I .46070, .02420) ,( 1 .43231, . 11618) , 

2( 1 .38757, .19766 ) , ( 1 .32928, .26620 ) , (1.26056, .31997), 

2( 1. 1846 v8,* 35775 ), ( 1.10497, .37397 ), (1 .02473 ,. 38 36 7) , 

2 ( .94712, .37250) ,( .87502, .34665) , ( .81 101, .30780) , 

2( .75 721 , .25 8 02) , C .7 1531 ,. 1 9972 ),(. 66647 1355 1 ) , 

2 ( .67131, .06814) ,( .66993, .00036) ,( .68187, -.06517) , 

2( .70622,-.12595) , ( .74160 ,-.17976 ) , (. 78628 ,-.22471 ) , 

2( .33821 ,-.25931 ) , (. 8951? ,-.28253 ) , ( .95464, -.29377 ) , 

2( 1.0 143 5, -.29295) , ( 1 . 07 1 9 , 28043 ) ,( 1 . 12509 2570? ) , 

2( 1.1 71 9 3, -.2239 3) , ( 1 .2 10 75 , - . 16270 ) , (1 . 2402 1 , 1 351 6 ) , 

2 { 1.25939,^,08335) , ( 1 .26778 ,-.02940 ) , ( 1 .26529 , ,02451 ) , 

2 ( 1.2 522 7, .0 7626) , (1.22947,0.1238 5) , ( 1 . 19799 , . 1 65 50) , 

2(1. 1592 6,. 19969), (1.1 14997,. 22523), (1.06701,. 241 29)/ 

DATA PI/3.14159/ 

data a/. 06233, .00404, .00101 , .00054, 

2.00040 . . 00028 , *00013 , .00003/ 

DATA S/.79786,. 012.56, .00179, .00067, 

2.00041. . 00025. .00011. .00002/ 

1 F( X.LE ..9) GO TO 200 . 

IF(X.L7.10.) GO TO 100 
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VJ=:.D 

VY=.0 

DO 10 1=^1, 6 
K = l-1 

VJ=VJ+A{ n=*(E*/X)**(2’>K-fl)*SlG*vj 
VY~VY+B( 1 )^(8./X)**<2 *K ) aSIGN 
10 CCA'Tlr^UF. 

VH=l •~C^^PlX(VJ,-VY)♦CEXP(CMPLX^ ,0^P1 /4*-X) )/SORT<X) 

RETURN 

100 CONTINUE 
y=5.*x 
J=Y+1,5 
IF(J.LT,2) 

IF(J,GT*50I J=EO 

JC>=J + 1 
FJ = J 

YJ=FJ-1 • 

Q=Y-YJ 

CT=0/2, 

C=OTy(Q~l.) 

D=QT*(Q^1,) 

7H=C*0( JMI+D*G( JP) +F^G(J) 

VH=CONJG(VH) 

^.TTUPN 

200 CONTINUE 

CAUL HANK{X,H0,Hlr2) 

HBO = X~X *^3/9. +X++5/ (9.^25* *^99. ) 

HBl =?Xyx/3.~X<'«4/4 5.+X’>«6/ (63.^25. ) -X** 8/ ( 8 1 .*25 • » 
VH=X«HO+X4(HBO«H1-HB1*HO) 

RETURN 

END 

SU3R0UT INt: CFF ( X A , YA , XB , Y5 »DK , CPH, SPH r ZS ? E JjEMjE JE) 
CGMPLFX ZStEJtE JA,FJ0,F,EM»FJEf F1,F2 
C0MM0N/C0A/CCC,EP2f TSK2L 
data PT/3, 14X59/ 

ETA=120*’^P1 
CA- (XB-XA)/0K 
CB=:( Y6-YA)/DK 
A=XA+CPH+YA*SPH 

b=xr*cph+yb«sph 

E JA=CMPLX (C09( A) t S1N( A) ) 

^^JB-CMPLX(CnS (8 ) tSlNIB) ) 

C=CA*CPH+CB^SPH 

S=Ce#CPH-CAySPH 

F=CMPLX{ ,0,DK )-^EJA 

I F< ABS(C) -GT<>eOOn ( E JB-F J A ) /C 

FJ=-30.«( ,707,".707)*F 

EM=2S« ( *707,-*707) (4,*P 1 ) 
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cc 


SP^l. + S/SORtf PR2) 

SM=i*-s/s{;'RT<rP2) 

Fl = C^»PLX( Oc,7SK2L ) 
f^^=CMPLX(0.t“TSK2L) 

IF< S.P.GT. 0.001 ) F2=(CFXP(CMPLX( .0»TSK?L«SP) )~1 .) /SP 
IF( SM.GT. 0.001) Pl=(CPXP(CMPLXr .0»-TSK2L*SM) )-l. )/SM 
e JE = *“F J «( r-R 2- i .)/ FR2#( F 1 + F2 5 

RETURN* 

END 

SUBROUT jNF CELSfXl ,Y1 rX2,Y2rXS,YS»DKt INTtPl) 

COMPLEX hQtHItPI 
DATA Pr/3.141F9/ 

CBET=(X2-X1 ) /OK 

SRET = < Y2-YD/DK 

XA=( XS-Xl )*CPFT+( YS~Y1 )’>SBET 

YA=-(XS-X1)=<^SBET+(YS-Y1 )’i^CBET 

X=XA 

Y=YA 

Y$O^Y*«? 

P.MIN=ABS(Y) 

lE(X.LT.O.O) PMIN=SORTCX*X-*-YSO) 

IF( X.GT .OK) PMIN = SQRT( (X~0K )^+2+YS0) 

ENT=1+(4*INT ) /I 0 
I S5=ENT«0k/RM 1N 
I SS=2*( 1 SS/2) 
iFf ISS.LT.2) 1SS=2 
F1T=ISS 
IS0=ISS+1 
OS=DK/FIT 
SGI=-1. 

XP=0. 

Pl=(0.t0.) 

DO 100 1=1,ISQ 
C=SG1>3. 

IF( 1 .FO.I.OR.r.FO.ISO) C=l. 

OELX=X~XP 

RK= SQRT(0ELX**2+YS0) 

CALL HANKCRKrHO.HlTO) 

P1=P1+H0«C 
SGi=-scr 
XP=XP+DS 
100 CONTINUE 

Pl = Pl<f r~30.*PI) «D5/3. 

RETURN 

END 

SUBROUTINE R E C ( A t B , 1 DM, NM, NP t X , Y , 1 A, 1 8 > 

DIMENSION X( IDM) ,Y( IDM) , I A( IDM) , 1B( IDM) 

NX=A*5.'^1.5 
NY=B*5.+1 *5 
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IF(NX.LT*4)^X=4 

IF(NY.IT.4)NY=4 

dx=a/mx 

DY=B/NY 
NM=?=^ (X»X + N'Y ) 

DO 3 J = 

1A( J)=J 
IB( J)^J+1 

IF (J .FQ. m) 1R(J) = 1 

NP = NK 
\!YP-NY+1 
1 T-2♦^!V-►^JX + 2 
AT=,5«A 
BT=.5*F 

no 1 1=^1 tNYP 
1 I- IT- I 
XU )-AT 

Y(I )=-RT-M I-l )+OY 
XU! )--X( 1 } 

Y(ll)-Y(n 
1 T=MYP + 2« ( ) 

IS=NY+.? 

IF=MX+NY 

no ? i = isuE 

I I-IS) 

X ( I ) = AT-n-I S + 1 ) + DX 
YU )=BT 

xun=x( n 

Y(in=“Y{j) 

RFTUFN 

END 



APPENDIX C 

COMPUTER PROGRAMS FOR RADIATION AND SCATTERING 
FROM TE DIELECTRIC-COATED CYLINDERS 


COMPLEX COT,ZS»Ill»Z12tYII fVl2»HZM»HZS*HZT 

COMPLEX C(56,56),CJ(56) ,H J J ( 56 ) tHMMi 56 ) ,V J < 56) tHJ JEQ( 56) 

niME\'SION lA(56>,I8(56),U(36>,12(56)f l3C56)f JA(56) 

OrMEK'SIOfv MO(56,5) ,NO(56) »X( 56)tY(56) , 0 ( 56 ) * XC ( 56 1 ♦ YC (56 ) ,DC(5f 
COMMON/CO A /TSK,ER2, CONST 
data IDMr I NT/ 56 f 10/ 

DATA PI tTP,ETA/3.14159,6.2«31Rf376,727/ 

7 FDRMAT( 1X*BE15.7) 

5 EOPMATIIHO) 

7 F0RMAT(7F10.5) 

8 FORMAT! IX, 1415) 

C0T=1,414214«ETA«CMPLX( l.,~) •) 

10 REA0<5,8) IWP,LOP*NM,NP 

WRITE! 6,8)1 WR, LOP ,NM,NP 
WPITF16,5) 

DO 50 J=1,NM 
READ(5, 8) IA( J) , IB(J) 

50 WR!TE(6,B)J,IA( J) ,IB( J) 

WR1TF(6,5) 

CALL SORT ( I A, i B , 1 1 , 1 2 ♦ I 3 , JA , JB , MO, NO , MM ,NP , N, 1 DM ,MA X t Ml N ) 
IF(M1N,LT.1 .OR. MAX.GT*5)Gn TO 300 
DC 35 1 = 1, N 

55 WRITE(6, 8)1, 11(11,12(1), 13(1) 

WR1TE(6,5) 

OD 60 1=1, NP 
RFAD(5 ,7)XC( n ,YC( 1 ) 

FI=I 

60 WR17E(6,2)FI,XCn),YCU) 

WRITE (6, 5) 

00 70 JAN=1,NM 
K1M=1A(JAN) 

L1M=IB(JAN) 

70 DC( JAN) = SORT ( ( XC ( L I M ) -XC ( K I M ) ) ( YC ( L I M )“YC ( K IM )) ’^★Z) 

READ (5, 7) rSL,r:R2 
WRlTn6,2> TSL,FR2 
WRITE(6,5) 

TSK=TP*TSL 

CDNST=( FR2~X. )/FR2*( COS (TSK*S0RT(ER2~1. ))•!.) 

80 RFAD(5,7)CMM,DPH,FMC,SCALE,TC 
WRI TE(6,2)CMM,DPH,FMCtSCALF,TC 
WRITE (6, 5) 

WAVM=300./FMC 

TPL=TP 



IF{ SCALP «GT.O. )TPL = TP*SCALe/WAVM 
no PO IAN'=1,NP 
K(1 AN)=TPL=3‘XC( 1AM) 

90 Y n AN) =TPL^YC ( IAN) 
on 95 JAN=^ 1 ,NM 
95 0( JAN)“TPL*nC ( JAN) 

112 = 1 

1SYH=0 

ZS=(,0,.0) 

tk=tpl*tc 

I F( CMM.GT, 0. )CALL C SURF ( CMM , FMC , TX , Z S ) 

IF(CMM.GT.O^) ISYM=1 
TSL.GT.O.O) 1SYM=1 

CALL C0ANT(C,D,X,Y»Z^,lA,lEtU,I2,I3,l$YM 
Bt I dm » 1 N T , JA , JB , MD, N , ND, NM , NP ) 

IF( 1SYM*F0. 10)GG TO 300 
GO TO Ul0fl20»130,l40),L0P 
110 9FA0(5,8)rGN 

CALL VAUSnOM,lGNtlSYM, lWRtI12»N»CtCJtYll ) 

FGN=1GN 

WRIlF(6t2)FGN,Vll 

GO TO 200 

120 RHAD (5,P) JSAtJSB 

CALL VWASdA, IB, IDM, ISYM, IWP , 1 1 , I 2, I 3, 112 , JSA, JSB ,MD, N,ND ,NM 
2,C,CJ,D,VJ,Y11) 

FSA=JSA 

FSB=JSS 

WRrTF(6,2)FSA,FSB,Yll 
GO TO 200 

130 RFAn(5,7)PSJ tXCSfYCS 
XS-TPL’^'XCS 
YS=TPl«YCS 

CALL VMLS( I A, IB, lOM, INT» 1 SYM , I WR , 11,12,13, 1 12 ,MD ,N, NO ,NI^ , 
2C,CJ,D,PSl,VJ9X,V,XSfYS,Yll,2S) 

WRITF(6 ,2)PSr rXCS,YCS,YU 

GO TO 200 

lAO PFAn(5,7)BSC,PMl 
WRITE(6,2)BSC,PHI 
200 WRITE (6, 5) 

in LOP»NC*A)G=ReAL(Yll) 

1NC=-1 

IF( L0P,EQ,4) INC=1 
1PA = 2 

IP( L0P*6Q.A) IPA=1 
IF( L0P.NF*4)BSC=-1. 

NPH=3C0*/DPH>1 .5 
GR=,0 

DO 280 IPH=IPA,NPH 

EPH=lPH-2 

PH=DPH’J^FPH 
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1 F( IPH.Pg.l JPHspHI 

CALL VFFHA, 18» I\Xt IDMtlSYMflWR, II ,12, 1 3, 1 1 2 , LOP ,M0 ,N ,MD ,NW{, 
2C tC J, D, t: WL, G, GA IN, HJ JtHMM, HZ StMZT, PH, ECS,VJ,X,Y,XS,YS,ZS) 

I F{ LnP.NE,4)WRlTE(6,2)PH,CAlM 
IF( LOP .EQ *4) WRITE (6, 2) PH, FWL,FCS 
I!XC=-1 

]F(BSC«GT.O.) J?^C = 1 
2SO G'?=GR+CA8S (HZT)*^2 
WRlTr<6,‘5) 

$CS=.0174533*DPH*GR 

IF( L0P.E0.4)WR ITE(6,2)SCS 

GR=FTA*SCS ‘ 

IF( LCJP.NE.4>W».ITE(6,2)GR 
WRITT{6,5) 

READ{S,B) jOBtLOP 
1 F( JOB.EQ*IO)GO TG 10 
1F( JOB. TO. BO) GO TO BO 
I F( J0B.FQ.300 )G0 TO 300 
GO T0<I10,I20,130,140),L0P 
300 CONTINUE 
CALL FXIT 
E NO 

GUeROUTlNF SORT ( 1 A , I B , 1 1 , 1 2 , I 3 , J A , JB ,M0, ND ,NM ,NP , N, MAX, M I N 
2,ICJ,INM) 

DIMENSION IAm,IBa),NO(n,MO(INM,4),JSP(20) 

DIMENSION nm ,i2(i),i3(n,jAa),jBm 

9 FORMAT! 3X , »MAX = * , I 5 , 3X , • M I N = •,I5,3X,*N = *f\5) 

1=0 

DO 24 K = 1,NP 
NJK = 0 

00 20 J=i,NM 

1 ND ( I A ( J ) -K ) ! I B ( J ) --K ) 

IF! 1ND.NF.O)GO TO 20 
NJK==NJK + l 
JSP(NJK)=J 

20 CONTINUE 
MDD=NJK“1 

IF!M0D.LE.C)G0 TO 24 
DO 22 IMD=l,MnO 
1=1 + 1 

IF! I .GT.ICJIGD TO 22 
IPD=IM0+1 
JAI=JSP(IMD) 
j A ! n = j A I 

JBI=JSP(IPD) 

J9! 1 )=JBI 
n!n = iAUAi) 

IB! lA! JAl).FQ.K)Iin) = ]B!JAn 

12(I)=K 

l3n»=U(JBI) 

IF! lA! JBI) .EO.K) 13! n = lB! JBi ) 

22 CONTINUF 
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24 CONTINUF 
M=T 

DO ?0 J=l,NKi 
ND( J )^0 
DD 30 K“l.,4 
30 MD(JtK]==0 
III=N 

I Ft N.G7 .iC J) I U^ICJ 
DO 40 1=U 1 I I 
J=JA(I) 

DO 36 L=li? 

NDt = 

K=1 

M=0 

32 NJK=MD(J*K) 

IF(MJK.NE*0)G0 TO 34 
M = 1 

MDtJ?K)=I 
34 K=K+1 

lF(K,GTo4>GD TO 36 
lFtM«E0»0)GD TO 32 

38 j=jfttn 

40 COMTINUE 
MIN=100 
MAX==0 

DO 46 J=1,NM 
NDJ = f\[D(J) 

lF(NDJ.GT,MAX)^tAX=NDJ ' 

46 IF(ISIDJ*LT*MIM)MIN=!^0J 

IFtMAX»GT.4 .OR. MIN.LT.l .OR. N.GT. IC J ) WR I TF t 6 ^9 ) M AX 1 M 

RFTURN 

END 

SUBROUTINE CD AN T ( C ^ Or X ^ Y 7 Z S v I A * I B , 1 1 ^ 12, 13? ISYH 
2, IDM?INT? JA? JB?MD?N,ND 7 NM?NP) 

COMPLEX ZS,Pll*P12»P?l,P?2,01l,QI2fQ21 ,Q?2?P( 2??) ?Qt2?2I 
COMPLFX CtlOMfIDM) 

DTMFNSiON X(1 ) ?Yn ) ,Dtn ?IA< 1) flBti) f JAtX) f JBt 1) 

DIMENSION IK 1 ) 1 12(1) 713(1 )?KDf IDM?4) ?NDt n 
2 FORMATt 3X, “DMAX » ? ElO .3 ? 3X ? ’OMIN = SE10.3) 

DO 20 1=1, N 
DO 20 J=1,N 

20 cn,j) = ( .o,.o) 

DMAX=.0 

DMIN=100. 

DO 7b J=1,NM I 

K=IA(J) 

L=IB(J) 

D( J )=SORT( (X( L)~X(K ) Y( L)-Y(K) J ♦*2) 

IFtOt J) -GT.DMAX)DMAX=n( J) 

25 IFtOt J) ,LT.DMIN)DMIN=D(J) 

DRAT=DM IN/DMAX 

IFfDMAX.LT. 3. .AND. DRAT .GT . .01) GO TO 30 
N=0 

WRITF(672)0MAXtDMIN 
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30 DO ?00 K = 1 
MDK = M0< K) 

KA=IA(K) 

KB=1B(K) 

DK=DCK) 
on 200 L=1,NM 
NDL=ND( L) 

LA=IA(L) 

iB=ie(L> 

DL=0(L) 

MIL = 0 

DO 200 II=1tNDK 

r=MD(K,in 

FI=1« 

IF(KR 12( 1) )GD TU 36 
IF<KP.F0.1UI))FI=-1. 

IS=1 

GO TO ^0 
36 IF{KA.ED.13< 

IS=2 

AO DO 2C0 JJ^l ,NDL 

J^MDdrJJ) 

I Ft 1 SYM.NE.OGD TO 42 
IF( I .GT.J)GO TO 200 
42 FJ=1. 

IF(LB.TC.I2(J))GD TO 46 
IF( LB,EQ. 1 1 1 J ) )FJ=*1 . 

JS=1 

GO TO 5 0 

46 IF(LA.E0.I3( J> 

JS=2 

50 I F(N1L*NF,0)GG TO 168 

IF(k*EO.L)GO TO 120 
I MDx(lA~KA)’>( L6-KA)’J‘( LA-KB J’^^tLB-KB) 

1 Ft INO.FQ.OIGO TO 80 
C SEGMENTS K AND L SHARE NO PDJNTS 

CALL ZMM3(X(KA) ,Y(KA) tXtKB) ,Y(KB),y(LA) tYtLA),y(LB) ,Y tLB) 
2DK,DL» INT,P(l,lUP(lt2) ,P(2, 1 ),P(2,2n 

GO TO 166 

C SEGMENTS K AND L SHARE ONE POINT (THEY INTERSECT) 

80 K.G-3 

JM=KB 
JC=«<A 
KF=-1 

INO=(KB-LA)«t KB-LB ) 

I Ft IND.NE,0)G0 TO 82 
. jr.=KB 
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KF=1 

JM=KA 

KG=0 

82 LG=3 
JP=LA 
LF=-1 

IF( La*FQ.JC)GO TO 83 

JP=LB - 

LF=1 

LG=0 

83 SGM=KF«LF 

CALL 2^*M2^X^JM),Y(JM),X( JC)»Y(JC),X(JP),Y( JP),ZS,DK,DLt 
21 0(1,2), Q (2, n, 0(2,2)) 

DO 98 KK=1,2 
KP= IASS (KK->KG) 

DO 08 LL=1,2 
LP=IABS(LL^LG) 

98 P(KP,LP)=SGN*0(KK,LL) 

GO TO 188 

C K=L (SFLF PFACTIDN) OF SFGMFNT K) 

120 CALL ZMMl(DK,2S,P(l,U,P(l,2)) 

o(2,2)=PU,n 

168 C(I , J)=C ( 1 r J)+F!«FJ=«'P( IS, JS) 

200 CDNTIMUt 
PFTURV 
FNJD 

SUPRQUTINE ZMM1(DK,ZS,P11,P12) 

COMPLFX ZS,H0,H1,PI1,P12 
C 

C0M«0N/C0A/TSK,FR2?C0NST 

C 

DATA PI/3*1A159/ 

CDK=COS(DK) 

S0K=S1N(DK) 

CALL HAW(DK,H0,H1,2) 

SDKS-S0K^^2 

CDKS=CDK**2 

PH =~2.*Hl>^CDK+H0’J'SDK + 2,^(.0,l.)=4^n. + CnKS)/Pl/DK 
Pl2 =-H0*CDK*SDK + H1=«^( l«+CDKS)-4««(.0,l.)=^CDK/PI/nK 
P11 = 15**DK^PU/SDKS^»( 1* + C0NST) 

P12 = 15,^OK^Pl2/SDKS=fi‘( l. + CONST) 

C 

T0K=2.*DK 

CTDK=rCCS(TDK) 

STDK=SIN(TDK) 

CCT=:SiNn SK«SQRT( ER2-1. ) ) /SQRT( FRP'l . ) 

Pll = PlX-^CCT*=(eR2“l#)/FR2*(.0,7.5)^(T0K^'ST0K)/S0KS 
PX2 = Pl2*-iCCT*( FR2-1. )/FR2*( .0,7.5) >* 

2(CDKJ<^(TDK + SrDK)+SDK*(I.-CTDK) ) /SDKS 
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O j Cl a 


c 


RS=CABS(2S> 

IF(RS*Le*0.)G0 TO 100 
CST=16,*P1*S0KS 
TDK =2. ♦OK 
7DK=C0S(TDK) 

TDK-SIN<TDK) 

n = pii‘^zs* (tdk-stdk )/csT 

12 =P12 -»-ZS^n l.-CTDK) ♦SDK+(STDK-TDK)^CDK)/CST 
100 RFTURW 
F ND 

SUBROUriNE ZMM2(XltYl,X2,Y?,X3,Y3,ZS 
2,DK1,0K2, INTtOll tC!?,Q?l tO?2) 

COMPLEX H0,Hi,HH0,HHI,SH0,SHl,0ll,Ql2,C21»Q22 
COMPLEX DHH0,DHH1 »OHO,nHl ,DSH0,DSH1 

COMPLEX SU,S12,S21,S22,Tll,Tl2,T21tT22,Yll,Y12,Y21 tY22 
COMPLEX DTn,0712fDT21,DT22,DYllfnY12tDY21*DY22 
COMPLEX ZS,RKH1,SX1,SX2,CCP,FUM,C0T 

COMPLEX CX1,CX2»DS11 fDS12tOS21,OS22 
COMPLEX DP11,DP12,DP21 , DP22 , P 1 1 , PI 2 f P2 1 , P22 
CDMMON/COA/TSK ,r:P2fC0MST 

DATA CCP»PI/(.0*. 63662) *3.141^9/ 

SDK1=:SIN(DKU 
SnK2=SIN(DK2 ) 

CDK 1=CDS (DKl) 

C0K2=CnS(0K2) 

CBET=(X2-X1)/DK1 
SeET=(Y2~Yl) /DKl 
XB=(X3-X1)^C3ET^( Y3-Y1)*$BET 
VB=-(X3-Xn^S0ET+( Y3-Yn^CBET 
CAL=<XB~DK1)/DK2 
SAL = ABS(YFi/DK2) 

CALL HAMK (DK2,HH0fHHl ,2) 

DHHO=DK2’>HHO 

DHH1=0K2^HH1 

C1S2=CDK.1«SDK? 

C1C2 = C0K.1^C0K2 
IF(CAL, LT,Oe )G0 TO 20 
1 F( SAL.GT..DA)GO TO 20 
CNT=~1S .♦CAL/S0K1/SDK2 

Call hank(dki ,ho,hi,2) 

DHO=DK1*HO 

DKS=0K1+0K2 

call HANK(DKStSH0,SHl,2) 

DSHO=DKS^SHO 

0SH1=DKS«SK1 

011 = CNIT*(CDK1«OSH1-C1S2<‘DHO-C1C2«DHI-DHH1 + CCP^COK2) 
012=CNT^<CDK2^DHH1-SOK2*DHHO-CCP+CDK1*OH1+C1S2+DSHO-C 1C2^DSH1 ) 
Q21=CNT«( SDK2^DHC— DSHl+COKZ^DHl+COKl^OHHl-CCP^'CICZ) 

Q22 = CNT^(CXS2^0HH0-C1C2^DHH1+CCP*C0K1“0HI-S0K2*DSH0-»-CDK2*0SH1 ) 
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on=oii*a.^-coNST) 

QI? = 012«( l.+CONfST) 

oaiso^i’j^d^+coNST) 

Q22=022’»tl •■♦-CONST) 

UPTURN 

S11=-DHH1-^CCP*CDK2 
S l?=-'SnK2=^OHHO-^CDK2*DHHl-CCP 
S21=(DHH1-CCP^C0K2)*CDK1 
S22~(SDK2*DHH0-CDK2*DHH1-^CCP)*CDK1 

DSll=nHH0-^CCP«SDK2 

DS12=^CDK2’^'0HH0-SDK2*DHHl 

DS21^=CDKl*(-OSn) 

DS2?=CDK1#1-DS12) 

DKS1=DK1*^2 

AL=ATAN2(SAL,CAL) 

PMIN^OKl 

1F(CAL*GP*0*)G0 TO 30 

RMIN=DK1*SAL 

DCR=-DK1*CAL 

IF( DK2•LT.DCR)RMIN=SQRTCDKS1■|■2•^DK1♦DK2«CAL■^DK2«^K2) 
FNT=1-M4>^^INT)/10 
I NP=FNT#DK2/RI^1N 
1NP=2*{ INP/2) 

IF(1NP.LT.2)INP=2 

IP-INP+1 

DT=DK2/INP 

TK=*0 

SXl=(.Ot.O) 

SX2=(.0,*0) 

CX1=(*0,.0) 

CX2=(.0,.0) 

SGI=-1* 

DO 90 1=1, IP 
0=SGI+3. 

IFd.tO.l .OR. I.FQ.IP)n = l. 

TKS=TK^TK 

RK=S0RT(DKS142«*0K1^TK«CAL+TKS) 

CALL HANK(RK,H0,HI,0) 

S1=SIN(DK2-“TK) 

S2=SIN(TK) 

SXl=SXlfSi>HO»D 
SX2 = SX2-^S2=*H0«D 

C1=-C0S IDK2-TK) 

C2=C0S(TK| 

CXl=CXl■^Cl♦HO♦D 

CX2=CX2^^C2>H0^D 



SGT=-SG1 
TK=TK+DT 
SXl=^SXl*DT/3. 
SX?=$X2*DT/3 . 
S21=$21-SX1 
S22 = S22^*^X2 
S12=51?+CDK1*SX2 
S11 = SH ^COKl’frSXl 

CXl=CXl«OT/3 . 

CX2 = CX2’^DT/3r 
OSll-OSn + CDKl*CXl 
OZl2-OSl2-^Cmi^CX2 
OS21=DS21-CXl 
nS22=DS22-CX2 


lH9-2^{ INT/2) 
ip=iNp-n 
j p= I p 

Tll=t(.Ot.O) 

T12=( .Ot.O) 

T21=(rO,*0) 

T22=(,0,»0) 

COT = ( .Of .0) 

RS=CABS(ZS) 

y ii = { .Of .0) 

Y12=(.Of.O) 

Y21=(.0f.G) 

Y22={.0,.0) 

Pll=(.Of.O) 

P12=( .0, .0) 

P?l = ( .0, .0) 

P22=<.0f.0) 

B = .0 

1F(AL.LT. .05)00 TO 210 
ALT=AL/2. 

CALT=Cns;( ALT) 

SALT=S1N( ALT ) 
PCP=(DK1-»'DK2)«CALT 
RSP = (0K2«DK1 ) >^SALT 
PMC = AT A N2( PSP, POP) 

SGI=“1. 

ph=-alt 
DPH=AL/IMP 
00 200 I- I, IP 
0=SGI+3. 

IF(I,FQ.l .OP. I. EQ. IP) 0:^1 
S AP = SI NM ALT+PH) 

SAM=SIM( ALT-PH) 



IFfPH.LE.PHC )RMAX=DK1>SAL/SAM 

I F( PH. GT*PHC )PMAX=^DK2»S AL/S A P 

DPX=:HMAX/1NP 

PK:r . 0 

SGJ--1 • 

0T11=(,0,.0) 

0T12=(.O,.0) 

DT21=(*D,.0) 

DT22=f •0,.0) 

DYl 1={ ,0, .0) 

OYl 2=( ^0,^0) 

DY21=( .0? .0) 

DY??=(*0,.0) 

DPI 1=( *0^.0) 

DP12=( .0,.0) 
nP2l=^( *0,*0) 

DP22 = { .0* .0) 


P.KH1=CCP 

no ino j=i,jp 

C ^SG J+3 e 

rF(J.FQ*l .OP. J.EQ.JP)C-1. 

1F(J.FQ,1)G0 TO 94 

CALL HAA’K(RK»HO»Hlt 1) 

RK,H1=RK*HI 

COIViT IMUE 

SK=RK* SAM/SAL 

TK=RK^SAP/SAL 

CX=COS(SK) 

C2=CGS tDKl-SK ) 

Sl^ SIN(aK2-TK ) 

S2=SIM(TK) 

FUN=^C*RKH1 
DYT l^DYl l-FUN^Cl«SX 
DY1?=DY 12-F/|M5?C1*S2 
0Y21 = DY2l-»-FU?^^C2*Sl 
Oy2 2 = QY22‘^FUN^C2’«'S2 


SlP=^-CnS(DK2~TK.) ^ 
S2P=C0S(TK) 

DPlls^DPli-FUM’fSlP’^Cl 
DP12=rDP12*-FUN^S2P'i»Cl 
0P21=0P2l + FUK‘*SlP^C2 
DP?2=DP22'^FUN!yS2P*C2 


SGJ=-SG J 
RK=RK-»-nRK 

1 F( RS.LF.O. )G0 TO 100 
SS1”SIN(SK) 

S S2=S IM (DK1“SK) 



om=oTn^FUN=»ssi*si 
0Ti;?=0Tl?+FUM’^SSl^S2 
DT 21 =DT 21 +FUN*SS 2 «S 1 
DT 2 2 :=DT 22 + FUNI*SS 2 ^S 2 
100 CONTINUE 

B=^SAP*nRK^D 
Yil=Yll+B«nYll 
Y12:=Y12 + S4^DY12 
Y21=Y21+B’«^DY21 
Y22 = Y22+B=^DY22 

Pll=Pll 43 ^nPll 
P 12 =P 12 +B*DP 12 
P 21 =P 214 B’^ 0 P 21 
P 22 = P 22 + B=«=DP 22 

PH=PH 4 DPH 

SGI=-SGI 

1F(RS*LF .0»)GD TO 200 
Tll=Tll4B’0^DTll 
T12=T124B*DT12 
T21=T21+B*DT21 
T22=T22+B^DT22 
200 CONTINUE 
B=0PH/9. 

IFIRS.GT.O. )COT:s( . 0 , 1 . )*ZS« 0 PH/( 72 **PI«SDK 1 ^SDK 2 *SAU 
210 CNT=~ 1 'S./SOK 1 /SOK 2 

QU=CNT*(CAL*( 1 ♦4CONST)^S11 + B*Y11)+COT*T1 1 
012 = CNT=J‘ (CAL* (1.4CONST)*S124B^Y12>4COT*T12 
021=CNT*(CAL*( 1 . + CONSTM^S21-^8*Y21>+COT*T21 
g22=CNT«(CAL’<^( 1.4CONST)^S224B^Y22)+COT^T2 2 

CCT=SIN(TSK*SQRT(ER 2 - 1 . n/S 0 RT(ER 2 -I .J 
QH=Qll-CN f^CCT«( ER 2 -l^) /FR 2 <‘( 0 S 11 *SAL~B*P 11 *CAL/SAU 
012 - 012 -‘CNT>CCT*(FR 2 -I* )/ER 2 «IDS 12 *SAL~B’«^P 12 *CAL/SAL) 
Q 21=02 3 -CNT^CCT*( ER 2 - 1 . )/FR 2 *( 0 S 21 »SAL--B»P 21 *C AL/SAU 
0 22 =Q 22 ~CNT 4 ^CCT*( FR 2 - 1 *) /FR 2 ’«'( 0 S 22 *SAL-B*P 2 ?*CAL/SAL) 

RETURN 

FNO 

SUBROUTINE 2 MM 3 ( X 1 » Y 1 , X 2 , Y 2 t X 3 t Y 3 » X 4 t Y 4 » Z S » 

?DKl, 0 K 2 » INT,PlltP 12 ,P 21 »P 22 ) 

COMPLEX HHA,HHB,ZS»HZltH 22 »CQT,ETl,ET 2 ,HO,Hl 

COMPLEX Pll ,P 12 ,P 21 ,P? 2 ,Sli,S 12 ,S 21 ,S 22 ,Tll,T 12 f T 21 ,T 22 

COMPLEX EQl,E02,OXl,012t021fQ22 
COMMON /COA/TSK,ER 2 tCONST 

DIMENSION CCl( 21 )»SSl< 21 )»CC 2 ( 2 nfSS 2 C 21 ) 
data ETATPI/ 376 * 727 t 3 -l 4159 / 

SU = U0».0) 

S 12 =(. 0 ,, 0 ) 

S21=(.Ot.O) 
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c 


S22=( 0) 

TU = (.0,*0) 

T12=( ,0, ,0) 
r?i-{ , 0 , . 0 ) 

T22=f.0,»0) 

Qll=( .0 , .0) 

012=( .0, .0) 

021 = 1 . 0,.0) 

Q22=(c0,.0) 

CPFT=(X2-Xl)/DK1 
SBFT=«Y2-Y1)/DK1 
XA:= (X3-xn>^CBET+(Y3-YU^SBET 
XB=(X4-X1 )*CBET+( Y4~Y1)^SBET 
YA = ~tX3~Xn*SBFT+ IY3-Y1)*CBET 
YB=-(X4-X1)’«'SBET+ (V4-Y1 ) YC6E T 
CAL = UB-XA)/DK2 
SAL=(YB-YA) /DK2 
RMIIVJ = lOG00. 

X=XA 

Y=YA 

OX=DK2*CAL/4, 

0Y=sDK2*SAL/4. 

DO 40 J =1,5 
YS=.0 
R = A B S ( Y ) 

IF(R.GT.l ,F-15) YS=Y=^Y 
xs=.o 

XA3=ABS (X-DKl ) 
r F( XAB.GT.l .£“15) XS=XAB*XAB 
I F{X.LT.O.)R=SQRT{y*X+YS) 
IF(X.GT.DK1)R=SQRT(XS+YS ) 
IF(R.LT.RMIN!)RMIN = R 
X=X+DX 
40 Y=Y+DY 

FM=1 + (4^IK!T)/10 
1SS=FNT*DK1/RM1N 
ISS = 2=4t( ISS/2) 

IF{ ISS-Lt.2)ISS=2 

Ifi ISS.GT.20) 1SS=20 

FSS=1$S 

ISO=ISS+l 

DS=DK1/FSS 

TTT = F?MT>DK2/RMIN 

ITT=2*( ITT/2) 

IFI 1TT.UT.2) ITT=2 
TF( ITT.GT.20) tTT=20 
FTT=ITT 

ITQ=ITT41 

DT=DK2/FTT 

XP=.0 
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SGN=-1* 

RS=CABS(2S) 

JUMP=:0 

ASAL=ABS(SAL) 

IF(RS.L?r#0* •ANJD.ASAL.lt. *04) JUMP=1 
IF( JUMP.FO.DGO in 60 
DO 50 1=1, I SQ 
C=SGN+3. 

IFfl.eQ.l .OR, I.EQ.ISQ)C=1 . 

cri( n=c«cos(DKi-xp) 

SSI ( n-C*SIN(DKl-XP) 
cc2n )=c*cos( xp) 

SS2 { I >=C*SIN(XP) 

SGN=~SG^J 
50 XP=XP+DS 
60 DX=DT»CAL 
OY=DT*SAL 
X=XA 
Y=YA 
TK=.0 
SGJ=-1. 

CDK1=C0S(DK1 ) 

DO 200 J=1,ITQ 
D=$GJ+3. 

IF(J.PQ,1 ,0R. J.P0.ITQ)D=1. 

CTl-D^S IM(OK?-TK) 

CT2=0’^SIN(TK) 

CTPl=-D>i‘CDS(0K2“TK ) 

CTP2=D*CDS< TK ) 

XP=.0 

YS=.0 

YA3=ABS(Y) 

i F( YAB.GT.1.E~15)YS=:YAP=4^YAB 
FTl = r.OfO) 

ET2=(.0,.0) 

»“'Z1 = (.0,.0) 

HZ2=(.0,.0) 

RKA=S0RT1X=«‘X+VS) 

RKB = SQR T( (X-OKi )«^^2 + YS) 
sph=yab/rka^yab/rkb 

IF(SPH.LT..04 .or. JUMP.FQ.llGO TO 110 
DO 100 1=1, ISQ 
OFLX=ABS(X-XP ) 

DXS=.0 

lF(DELX.GT.,l.F-151DXS=Dr:LX*nELX 

RK=SORT(DXS+YS) 

5PH=Y/RK 

ci=ccim 

S1=SS1(I) 

C2=CC2(I) 

S?=SS2( I ) 
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CALL HANK(RK,HOtHl,l) 

Fn = ETl-*Cl’frSPH«Hl 
?='T?=reT2 + C2’!‘SPH’^Hl 
XP=XP+D£ 

IF( RS,LL*0, )GQ TD ICO 
HZ.l=HZl+Sl*Hl=«‘SPH 

ICO CnMTINUE 

110 CALL HANK(RKA,HHA7H1,0) 

CALL HAMK(RKBtH0,H1 ,0) 

C 

E01 = ET l ’>CAL<^DS/3»-SAL*( CDK 
F0R = ET?^CAL*DS/3.-SAL=^^{CDK1«H0-HHA» 
Qll=0il+CTP1*EQ1 
QI2:=Q12«>CTP?^FQ1 
021=Q21+CTP1*E02 
Q22=022+CTP2*HQ2 
C 

FTl=FTl«SAL’^DS/3.+CAL’^(C0Kl»HHA“H0) 

ET2^ET2*SAL + DS/3*+CAL=fl^(CDKl^H0-HHA) 

S11=S11+CT1^ET1 

S12=£l?+CT2trTl 

S2l=S21+CTl*ET2 

S22=S22+CT2*ET2 

CGJ=-SGJ 

TK=TK+DT 

X=X+OX 

y=Y-»^DY 

IF(R$.LE.O.)GO TO 200 
m = Tll + CTl’S‘HZl 
T12=T12+CT2*HZ1 
T21=:T21+CT1=^H22 
T22 = T224'CT2=«'HZ2 
.200 CO^^T^NUF 

S0K1=:SIM(DK11 

SDK2=SrN(DK2) 

CST = -“ETA^DT/( 24.=J'PI^SDK1’^SDK2) 

CQT = (,0,1,)»DS=^DT^ZS/(7 2,*PI^SDK1^SDK7) 
P11=CST*( X,+CON STIR'S 11 + COT*! 11 
P 12=CSr^( 1,+CDNST )+S12+CDT*T12 
P21=CST«( X.+CONST)*S21+COT^T21 
P22 = C£T4f{l, + CDNST KS22+CQT#T22 
C 

cc 

CCT = SI?^< TSK*SQPT (ER2“1. ) ) /SQRT< ER2“1 • » 

CCT=CCT’^( « )/ER2 

PI1~P11-CC1^CI1=«^CST 

P12 = P12-CCT>S‘012ACST 

P21=:P2l~CCT*0 21*CST 

P22^P2 2-CCT=^^Q22^CS7 
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RfTUP.W 

END 

*;UBROUTINE HANK(X tH,Hl, ID) 

COMPLEX H»H1 

DATA TSP/. 63661977/ 

IF(X*GT,3.)60 TO 100 
XtN=TSP»AL0G(X/2.) 

B = .0 
B1=.0 
Y=*0 
Y1=^*0 
Xl = X/3. 

X2=X1+X1 

IE(X1«LT,.1)G0 TD 60 

XA=X2^X2 

X6=X2>^XA 

IF(X1.LT*.3)G0 TO 55 

X8=X2’^^X6 

X10=X2*X8 

X12sX2«X10 

B = .21E-3*X12-,39444F--2«X10^*A44479E~1*X8 
Y=-.24846E-3«X12-*'*427916E“2*X10-.4261214E-1*X8 
B1=.1109E~4*X12-.31761E~3*X10'^.443319E-2«X8 
Yl= .2787?E-2>>^X12-.400976F-14X10+.31239 51»X8 
55 B=S-.3l63866’f^X6 + l*2656208^X4 

Y=Y+.25300117*X6-.743503S4«X4 
8 1=^8 395428 9E-l*X6’f. 2 1093 573^X4 
Y1=Y1-1.3164827*X6+2.1682709’»^X4 
60 B=R-2.2499997*X2+1, 

Y=Y+.60559366*X2+.36746691+XLN«B 

Bl=X*(Bl-.56249985*X2+.5) 

Yl=(Yl+*2212091*X2-. 6366198 »/X+XLN*Bl 
GO TO 2 0.0 
100 SW=SQRT<X) 

Xl=3./X 

X2=X1*X1 

X3=X1«X2 

X4=X1<^X3 

X5-X1+X4 

X6=X1*X5 

F=.797B8456-,77E-6*X1~,55274E-2’^X2-.9512E~4*X3+,137237E-2’^X4 

2~.72805F-*3*X5+.14476E~3*X6 

T=X~*78539ei6-,4l66397E-l^Xi«.3954E-4*X2+*262573E-2^X3 
2“*54125E-3*X4~*29333E-3*X5^*13558E-3»X6 
B^E^frCOS (T) /SW 
Y=F^SIN(T)/SW 

F =. 7978 8456 156F-5^XH“. 16 59667E-1*X2+. 17 105E-3>X3“. 24951 1E-2*X4 

2+,113653E-2*X5“.20032E«-3’<^X6 

T=X«-2*3 561945+.12499612=«'Xl + *565E-4*X2~*637879E-2»X3^*7434eE-3*X4 
2+ .79824E-3^X5“«29166E~3«X6 
B1=F+C0S(T)/SW 
Yl=F*SIN(T) /SW 
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200 M=CMPLX<B»-Y) 

M1=CMPLX { B1 t-Yl ) 

PFTURN 
E NO 

SUB ROU time CRCUK C»S* ICC , i SY M , IWR 1 1 1 2 » N ) 

complex cncc*icc)9su) 

COMPLEX F»P,SS,T 

2 FORMAT ( lX,XI5,lF10«3»lFlS.7aE10^0) 

5 FORMAT(IHO) 

IF< 112.ME.DGD TO 22 
I F( M«FQ* 1) S(I )=S( 1 ) /C( 1,1) 

IF(M*FQ*1)GD TO 100 
IF(TSYM.MF.O)GO TO P 
00 6 I=1,M 
DO 6 I 

6 c f j , 1 )=cn , j) 

fi F=:C(1,X) 

DO 10 L“2,M 

10 C (1 ,L) = C (1 ,u / F 

00 20 L=2,M 
LLL=L“1 

DO 20 1=L,M 
*'■ =C (I , L) 

DO n K = l,LlL 

11 F^F-C( I ,K)«C(Ka) 

C(I,U-F 

1 F( L »EQr I) CD TO 20 
P=C(L,L) 

IF( ISYM*EQ.O)GO TO 15 

F=C { L, n 

DO 12 K^ULLL 

12 f"=F“C(L,K)«r (K,I) 
ca,n=F/p 

GO TU 20 
15 F=rC(I,L) 

C(L,I)=F/P 
20 COMTIMUE 
2? on 30 L=1,M 
P:^C(L,L) 

T=S(L) 

IF(L.pq* 1)G0 70 30 
LLL=L“1 
on 25 K = ULLL 
25 T=T“Ca*K)*$(k) 

30 S(L) = T/P 

no 38 L=2,M 

1 1 = 1+1 

r=s < I ) 

00 35 K=I1 ,M 
35 T = T-C < I ,K)«S (X ) 
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38 S(n=T 

1F(IWR,LE.0) GO TO 100 

CN0R=.0 

00 40 

SA=CABS(S(n ) 

40 1 F( SA.GT.CNORICNOR=SA 

IF(CN0R*LE*0 . )CNOR:tK 
DO 44 1=1 

ss=sn ) 

SA=CABS(SS> 

SNOR=SA/CNOR 

PH=.0 

IF( SA.GT*0.)PH=57*29578*ATAN2( A1MAG(SS),REAL(SS) ) 

44 WRITF(6t?) If SNORtSA, PH 
WRITE(6f5) 

100 RETURN 
END 

SUBROUTINE VNAS( IDMf IGNf ISYMf I WR f 112 fNfCfCJfYll) 

COMPLEX CnDM,IDM),Cjn)fYn 
00 20 I=lfN 
20 CJ(1)=(.0*.0) 

CJUGNr = (l.,0.) 

CALL CROUKCfCJflOMf ISYM, iWRf il2fN) 

112=2 

Y11=CJ(1GN) 

RETURN 

END 

SUBROUTINE VWAS (I A, IBt IDM, I SYM* IWR » I If 12* 1 3 f 112 , JSA » J SB , MD* N ,ND 
2fNMf CtCJfDfVJfYll ) 

COMPLEX C(IDMflDM)fCJ(i>,VJU)fYll 

dimension I A( 1 ) f IB< 1) *11 ( 1 ) * i2( 1 ) f 13(1) *MD( lDMf4) *ND( l> fO ( 1 ) 
AK=.0 

DO 20 K=JSAfJSB 
20 AK=AK+0(K) 

DO 30 I=lfN 
30 VJ( n = ( .Of.O) 

IF( JSB.GT*JSA)GO TO 200 

K=JSA 

DK=D(K) 

V=( l.-COS(DK) )/( AK*SIN(0K) ) 

KA=IA(K ) 

KB=IB(K) 

NDK=NO(K) 

DO 140 11=1, NDK 
I=MD(K,li) 

FI = 1. 

IF(K8.E0.12( I ) )G0 TO 136 
IFiKB.EC.Il ( I ) )FI=~1, 

GO TO 140 

136 1F(KA.E0.I3( i) )FI=-lf 
140 VJ(I )=VJ(n+FI<‘V 
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GO TO 2 HO 
?00 KA=IA(JSA) 

Kn=lB(JSA) 

la=ia(j<:a>i) 

LB=1B(JSA^1) 

1N0=^(LA“K6 )«( LB-KB) 

I F( IMO* FC-«0)GO TO 210 
KA=1B(JSA) 

K B= 1 A ( J G A ) 

210 DO 2.B0 K = JSAvJSB 
0K=0(K) 

V=(1*“C0S<0K) )/(AK*SIN(DK)) 

NDK:^ND(K) 

00 ?40 n^l9NDK 

i=MD<K, in 

Fm. 

IF(KB.F0n2( 1 ) )G0 TO 2.36 
IF(K&*F0«U(n )F1=-U 
GO TO 240 

236 IFf KA.F0^13U ) ) F I=“i« 

240 V J( I )xvjn ) 1*V 

IF(K.eO.JSB)GD TO 250 
LA=IA{K+1) 

LB=I3(K4lJ 

KA=KB 
KF. = LA 

IF( LA* EO.KA)KB=LB 
250 C0MT1MU»^ 

280 DO 300 I=ltN 

300 cj( n=vj{i) 

CALL CROUT(C,CJ, IDM* ISYM,IWB, U2,N) 

112=2 

Yll=(.0,.0) 

DO 400 1=1, N 
400 Yll=Yli + VJ(I)*^CJm 

END 

SUBROUTINE VMLS C I A , I B , 1 DH , I NT t I S YM , I WR , 1 1 , I 2 , 1 3 ,112 , MD , ND , NM ^ 
2C ,C J,0,RSI ,VJ,X,YtXS,YS,Yll, 2S) 

COMPLEX C{ IDM, IDM) ,CJ(X ) , VJ( I ) tYll,Pl, P2,01,02,ZS 
DIMENSION lAdinBamUl), I2(1),I3(X> 
dimension MD(l0H,4),ND(n,Xm,Y(l),0(l} 

DATA ETA, TP/376. 727, 6. 28318/ 

DO 100 1=1 ,N 
VJ(1)=(.0,.0) 

100 CJ( I )=( .0, »0) 

DO 240 K=1,NM 
KA=IA(K) 

KP=IB(K) 

CALL CMLSIPSI ,X(KA» , Y(KA) ,X( KB) ,Y(KB > , XS,YS,D(K) , !NT, Pl,P2,Ql902 

D1=ZS*Q1 

02=ZS’9^Q2 

NOK=ND(K). 
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DO 240 I 1=1, NDK 
I=MD(K, in 
FI=1. 

IF(KB.FC. J2( I ) ) GO TO 236 
IF(KS,tQ.U( 1 ) ) FI=-i, 

cj{ n=cj( n+FnPi 
vj( T )=vj( T ) + F n(Pi4-on 

GO TO 240 

236 IF( KA.FO* !3{ n ) Fl=-1* 

CJ{ I )=C J( n+FJ*P2 
VJ( 1 >=VJ( i )+Fi*(P2+Q2) 

240 CONTI NiUt 

CALL CROUT(C,Csl,10M,lGYM, lWR,n?,N) 
fl2 = 2 

Yll^CMPLXdP/ (4.*ETA ) ,0. ) 
on 300 1=1, N 
300 YU=Yll + Cjn)*VJ(I) 

RETURN 

'^ND 

SU3R0UT INE CMLS(PS1 , X 1 , Y 1 , X2 , Y2 , XS , Y S, DK , I NT, PI , P 2, Ql ,02) 
COMPLEX CST,H0,H1,P1,P2,Q1,Q2 
DATA fcTA/376.727/ 

DKH = DK/25.. 

D1=$QRT ( (X$~XI)«’^2+(YS-Y1 )^*2 ) 

D2=SQRT ( (XS~X2 YS-Y2) 

P2=( .0, ,0) 

P1=CMPLX( .5*PSI/360.,0.) 

01=(»0,^.0) 

02=( .0, .0) 

IF(Dl.LT,DKH)Gn TO 200 
P1=(.0,.0) 

P2=CMPLX( .5^PSr/360.,0. ) 

IF<D2,LT*DKH)&0 to 200 
SDK = S1N(0K.) 

P1=(.0,.0) 

P2=(.0,.0) 

CBFT=(X2-X1) /DK 

S3«^T=( Y2-Y1) /OK 

XA= ( XS-X1)*CBET4( YS-Y1)«SBFT 

YA=-(XS~X1 l*SBET4<YS-Yn*CeFT 

X=XA 

Y=YA 

Y5;g=Y«^2 

RMIN=ASS(Y) 

I<=( X,LT 4 0. )RMIN = SQRT (X’^X^YSO) 

1F( X.GT.nK )RMIN=SORT( (X-DK)**2+YS0) 

FNT=14 (44^1NT ) /lO 
1 5S=FNT4DK/RMIN 
IF( 1SS.LT.2)1SS=2 
DS=OK/ISS 
XP=DS/2* 
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00 100 1 = 1 , iss 

0gLX=X“XP 

RK=S0RT (r»‘ELX^^2+YSQ) 

SPH=y/rk 

S1=SIN(DK-XP) 

S2=SIN(XP) 

CALL HANK (RK, HO, HI, 2) 

P1=P1 + S1*H1’>SPH 
P2=P2-^ S-2«H1*SPH 
01=01 + S l^^wo 
02=Q2-»-S2’J‘H0 
100 XP=XP+DS 

CST = ( *0, K ) *DS/|4.5$^SDK} 

P1=CST^P1 

P2=CST*P2 

CNT=~0S/(4e*ETA«^SnK> 

Q 1 = CNT*Q1 
Q2=CNT*02 
200 RFTURN 
FND 

SUB ROUT INF VFF ( I A , 1 8 , INC , J DM , I SYM, 1 WR, 1 1 , 1 2 , 13 , 1 1 2, LOP, MO , NO ,Nw 
2C, C J, 0, F WL, G, GA I N ,H J J, HMM, HZ S, H2T, PH, PCS, VJ,X,Y,XS,YS,ZS) 

COMPLEX C0T,D0T,HJ1 ,HJ2,HM1,HH2,HZM,H2S,H2T,ZS 
COMPLEX CJU ) ,HJJ( 1 ) ,HMM( 1) ,C( lOM, IDM) ,VJ ( 1) 

COMPLEX HJIJ,HJ2J,HJJFGK0) 

DIMENSION I AU ) , IB ( 1) , II (1) , 12(1 ) , I3C 1 ) ,MD( 

DIMENSION NDn),X(U,Y(l),D(l) 

COMPLEX CONPH 
COMMDN/COa/CONPH 
C0MM0!M/C0A/TSK,ER2,CnNST 
data FTA, TP/376. 727,6.2831a/ 

C0T=1.41421A*ETA^CMPLX< l.,~l,) 

I 1 SYM.NE «0) DOT=CQT*CONJG( ZS I /Z S 
FCS=.0 

PHR=r,o 174533* PH 
CPH=COS(PHR) 

SPH=SJ N( PHP ) 

DO 232 1=1, N 
HJJEQC i )=( »0, .0) 

HJJ( I) =f *0, .0) 

232 HMM( n = C •0,.0J 
00 2^0 K=1,NM 
KA=IA(K) 

KB=IB(K) 

CALL CFF (X (KA ) ,Y( KA ) , X{ KB) , Y (KB ) ,D (K ) 

2 , CDH,SPH,Z$,HJl,HJ2,HMi,HM2,HJlJ,HJ2J) 

NDK = N0(K ) 

00 260 1 1=1 ,NOK 

I=MD(K,in 

FI = 1. 
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IF(KB,FC*I2(1))G0 TO 236 
1F( KBdHQ.ll ( 1) )F 1=-1 . 

HJJ ( U=HJJ( n+Fl *HJ1 
HMM( J ) =HMM( 1 ) + FI’>^HMl 
HJJEO( n=HJjFQU 
GO TO 250 

236 IFf KA,F0*13fI ) >F1=~1* 

HJJ( n=HjJ( n+FI*HJ2 
MMM { 1) =HMM ( I )+F 
HJJEQ( n=HJjF(M IT-^F1*HJ2J 
250 CONTINUE 

I FT INC.LE.O )G0 TO 270 
DO 260 I=1tN 
CJ( I )=CQT«HJJ( 1 ) 

vj( n=cj{i ) 

260 1 F( ISYM«NE.O)Vjn )=VJ( I )-DQT*HMM{I ) 

CALL CBOUT(CtCJ, IDW, I SYM , 1 WB , I 1 2 ) 

112=^2 

DO 265 1 = 1, N 

265 ECS = eCS+BEAL(VJ(n*CONJG(CJC I))) 

FCS=ECS/ETA 
270 H2S=t*0,.O) 

DO 360 lsl,N 

360 HZS=HZS + CJ( I I^^THJJI I ) >-»*CJ( I ^’►(CONPH^HJ J( I I ) ) 

HA5=CABS(HZS) 

1F( LOP. FQ,4) FWL=TP*NAB»^HAB 
MZT=HZ5 

1F(L0P.NE.3)G0 TO 400 
PSI=XS>!^CPH+YS=frSPH 
HZ.'IsCMPLXlCOSfPSUtSINI PSl n 
H2M=-{ I .,1. )«HZM/(2.«1.414?145»FTA> 
m2T=HZS+H2M 

400 HAe=rCABS(H2T ) 

I F( L0P.LT.4) GA TN=TP^ETA*HAB*HAB/G 

PETURN 

*^N0 

SUBROUTINE CFf ( X A , Y A , XB , YB ,D K> , CPH, SPH, 2 S , HJ 1 ,H J2 ,HM 1 , HW2 , 

COF^PLEX EJA,EJB,CST,ZS,HJ1,HJ2,HM1,HM2 

COMPLEX El, F2,CSTT, PHASE, FBI, FB2,HJ1J,HJ2J 

complex conph 

common/chb/conph 

C OM HD N / C D A / T S K , E R 2 , CO NS T 

DATA FTA,PI/376.727,3.14I59/ 

C A= ( XB'-XA) /DK 
CS={YB-YA )/0K 
G=C A*CPH+CB4$Ph 
P=CB’!‘CPH->CA*SPH 
i;;K = p3^«2 

A=:XA*CPH + YA^SPH 
B^rXB’i^CPH+YB + SPH 
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ejA=CMPLX(COS( A) *S1N(A) ) 

FJB=CMPLX(COS(B) 7 SINIB) ) 

SnK=SIN(DK) 

CCiK = CGS(OK} 

IF(GK.LT..001 )G0 TO ?30 

CST = CMPLX|1 » ♦!*) /<4«.x^PI^SDK*l »414214^GK) 

HMii=CST4' { EJA*CMPLX(CDK,G^SDK)~EJB ) 

HN>2=CST’^ (FJB*CV!PLX(CDKv-G’»‘SDKJ-EJA) 

GH TO 300 

2 50 CSr = CMPLX(“l • a . ) /( ,414214^SDK) 

iF<GaT,oaGn to ?eo 

HM1=CST’MDK*E J 6 --SDK*ejA ) 

HM2 = CST=^‘( SDK<‘£ Je-DK-^^EJA ) 

GO TO 300 

280 HMl =CST« ( SDK^t JA-OK^EJB ) 

MM2=CST^ (DK^F JA-SDK*EJB ) 

300 HJ1=^P*HM1 

HMl=-2S+HMl/eTA 

HM2=“ZS*HM2/FTA 

C 

D=TSK/2./PI 
E1=CMPLX(.5«D,,0) 
e2=CMPLX(-.5«n, , 0 ) 

SP^ 2 .^PI’«‘(S 0 PT(FP 2 ~ia+P) 

SM=2.’«^P I*(SORT( FRa-'ia-P) 

SSP = SINUSP>J^D) 
c sp~cns( sp+D) 

$SM = SIN( 

CSM = COS( SM>«=D) 

IF( SP.GT,0«001)F2=( .Of~.5)=«'(CHPLX(CSP»-S$P?“l, )/SP 
IF( SW.GT.0*001 )FI = ( .Ot-.5>^MCMPLX(CSM9 SSM)~1 a /SM 
PHASF=E1~?'2 

CONPH=( *^0,ia«(ER2-ia/tR2*(-El--E2) 

IF(GK,GT.C,001 )G0 TO 500 
FBi = *5^(DK-FJA + SDK*E JB)><^PHA5E 
FB2= (DK’»'F JB + SOK’^e JA)*PHASF 

GO TO 400 
500 CONTINUE 

FR1 = PHASE^(CMPLX( .0 * G ) ••^^F, J 8 *«-CNPLX (SDK , ~G*C DK ) ’5^‘F JA ) 
F32=PHASE*(CMPLX( .OtG)^EJA-CMPLX(SDK ,G^CDK)^FJB) 
400 CONTINUE 

CSTJ=C:NPLX( 1 . , l • ) / (4.*PI*S0K«1 ,414214) 
1F(GK.GT,,001) 

?CSTJ=CMPLXn . ♦ ! . ) /(4.*PI>SDK^1 . 41421 4«GK) 

00 = 1 , 

HJ1J=00«CSTJ^( ER2-1* ) /FP 2<^G^ (-FBI ) 

HJ2J = 00«CSTJ*(ER2-I, ) /ER2’i=G* <-FB2 ) 
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RETURN 

END 

SUBROUTINE CSURFICMM, EMC,TK,Z$) 

COMPLEX ETA ,R,2StETST 

DATA Et FTA0,TP,U/P,8B433E~12f376.727t6.28318, 12.56^^4^ 
ALPH = SQRT(TP’5‘FMC*U*CMM/2 , )«1 .E6 
SOT^SQRTI TP=»=FMC^U/( > 

fTA=CMPLX(SQT,SQT) 

TAT = 2.<'TK=«'SQRT(CMM/(?.’{‘E>TP*FMC) ) 

2S=FTA 

IF(TAT.GT,61>.)G0 TO 100 
ETAT^EXP (“TAT ) 

ETBT=^C«PLX(C05( TAT) ,-SlN(TAf ) ) 

R“FTAT*ETBT^(ETA0“FTA )/ (ETAO+FTA ) 

ZS=ETA=J^(1.+R ) / ( 1»-R ) 

ICO RETURN 
END 
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APPENDIX D 

COMPUTER PROGRAMS FOR RADIATION AND SCATTERING 
FROM PERFECTLY-CONDUCTING PLATES 


OIMFMS IO\ C, (31) AO) y 40 
COMPLEX Z<AO,^D) ,ZntZ179/SS9ZIJ 
CQHPLFX VA»Ve?VC»VO»VTA t VTB , VPA , VPB » F S T 
COMPLEX VT(AO) 9VP(A0)»ET(A0) ,EP(A0) 

I DM=40 
P 1=3.14159? 

TP=2.^PI 

PFAD(5,9B)ZKL»TK,M$9.NMAXtNW 
96 F0RMAT(?FI0«,4,3I10) 

333 CONTINUE 
ZL=ZKL/PI 
TK=TK*TP 
NZ=NW 

CALL S1MWC(NMAX»C) 

NP=^NS-1 

N2S=NZ“1 

MOnFS=NW^NZS 

M0nFT-2=^M0DES 

AL=ZL<^?,/NZ 

HL=AL/NS 

DWL=ZL/\*Z 

DDW=DWL/ ( NMAX-1 ) 

00 11 K = lrNV^ 

DO 11 L=1,NZS 
Z 1 2 = ( 0 ® 0 , 0 * 0 ) 

N1=K-ML“1)*NW 
00 30 i l=ltl 
Y1=DDW=4‘( U-1)^TP 
Y7=Y1 
Y3=Y1 

no 30 JJ=1,NMAX 
DT=ATAN?(DWL, 1.414) /{ NMAX~1) 

CC=C(JJ) 

TH=.OT=>( JJ-1 ) 

DV=1 „41A^DWL/{NMAX~1 ) 

V-DV«( JJ~1) +(K-2)^0*707«DWL 
YA^V’^TP^l.AlA 

IF( K.EC.1.AND.L.LE.3) YA = TA!M { TH)*TP*1 .414 

ye=ya 

YC=YA 

2 11 = ( 0 . 0 , 0 . 0 ) 

00 10 l=l9NP 
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FI=C0S(i>I*2l/AL) 

21=(-ZL>0.5+( I-l) WL)*TP 
Z2=^Zl-»-HL«TP 
2 3=Z2+HL«TP 
DO 10 J=^1,NP 
2 J=~AL^O»5 + HL’«'J 
Fj=CnS( 2 /AL) 

ZA-C-ZL’^O.^+i J-1)*HL)^TP + ( L~1 ) ’>=2L’«^TP/MZ 

ZB=ZA+HL’^TP 

ZC=ZB+Ht ♦TP 

ZSS = ZI J(TK,TKtTK,YltY2,Y3,21t22,Z3,*0,.0,*0,YA,YBrYC, 2A,ZPtZC) 
ZYFi^Fj 
Z11 = ZU+ZS5 

10 CONTINUF 
Cl=K>0WL-1.^14^V 

I F( V*LT,(K-1 ) ♦0r707^DWL ) CL=1 . 414*V~ (K“2 ) *DWL 
lF{K.Nf ,1 ) ZlZ^^Zl^ + Zll^CC^DV^CL 
T F( K.EQ.l ,A\'D.L*GT,3) Z 12 = 212 + Z11*CC*DV*CL 

IF(K«^FQ.1*AND*L^LF^3) Z12=^Z12 + Zn*CC^DT*(DWL«-l .41A^TAN( TH)) / 
2(C0S( TH)*CDS( TH) ) 

30 CO^niNUF 
CK=2 . 

IF(K .FO.l .AND,L*LE.3) CK:^^. 

2 (1 tNI ) = Z12*CK.*0.707/3./(DWL*DWL> 

11 COMTINUF 

DO 40 M = 1 ,?^2S 

00 40 N‘=1-,NW 

1 =N-i- (M“U«NW 
MM{ 1 ) =M“1 
NN(1)=M-1 

40 CONTI^mF 

DO 1 1=1, MODES 
DO 1 J=1,M00FS 
^1 = 1 ABS (MM( in 

N2=l A6 5 (NN( J)-MN( 1) ) 

DO 2 M=1,MZS - 

DC 2 N=1,\'W 
N I=N+ (M-1 )^MW 

1F(NI1.FQ.M-'1 , AND.Nf2.FQ*M~l) GO TO 3 

2 COf^TINUE 

3 CONTINUE 

zn»j) = z(i,Ni) 

1 CONTINUE 

DO 33 I=1,MQDFS 
00 33 J=l tMODFS 
n=I+MODrrS 
JJ=J+MODES 
33 Z(II,JJ)=zn,J) 
c THIS PART FOR CROS S-COUPLING 


117 



Dnw=nWL/(NMAX“l) 
call SiMWC(MMAX,C) 

DO 2? K-1,NW 
DO 22 L-lf^fZS 
N)I=H>t L-1 )«MW 
DO 22 

DO 22 N=1,K!7S 
<NJ=M'»-{N-1 )=«'NiW-»'MGDFS 
Z 12=10.0,0.0) 

'10 31 Il = l,Nf^AX 
CI=C(I1) 

Xl-TK 

X2=TK 

X3=TK 

Yl=Onw^( Il~l)’^TP->-(K-'X)+DWL*TP 

Y2-Y1 

Y3=Y1 

DO 31 JJ=1*MMAX 
0 J=C ( JJ ) 

X A=0 .0 
XB-0.0 
XC=:0.0 

2A=TP«(7L“DDW*( JJ~1 )“(M-1)^DWL) 

2B=ZA 

ZC=ZA 

2 1 1 = { . o , 0 a 0 ) 

DC 32 I = 1,NP 
2 r = ~O.B=^AL + HL^M 
F I=cns( ZI*PI/AL ) 

71= I I-l )^HL=«‘TP+f L-D’ATP^ZL/NZ 
Z2 = Z l'»-HL*TP 


Z 3 = 2 ?^HL*TP 
no 32 J=1,NP 

z j=-o.b^al+hl«j 

P J=C02( ZJ+PI/AL ) 

YA= ( J-1 )’^HL*TP+(N-'1 )=«<TP’^ZL/NZ 

YR = YA*^HL’^'TP 

YC=YB+HL«TP 


7£S=ZIJ (Xl,X2,X3»Yi ,Y2,Y3»Z1 ,Z?,23,XA, XB,XC,YA,YB,YC, ZA, ZB,ZC ) 
2 «F 1 ’^PJ 
Z11=ZSS+Z11 
32 CGXaiNUF 

Z12 = 212+Z11’5'CIVCJ 
31 CONTl^iUF 

Z (^11 tNJ)=21?AnDW’«'nnw/(nWL*DWL) / ( 3.^*2) 

WP I TL ( 6 , 99 ) N I , Ni J , Z ( N 1 , N J ) 

99 FORMAT (2UG,2F10. 9) 

22 CQMTINUF 


118 



C THIS PART fl\a THE VOLTAGES 

CALL S 1MWC(NMAX,C) 

DDW=^DWL/(NMAX-1) 

THI=9C« 

PHI =0*0 

loo CO^^TIr^UE 

CTHI=COS (THI^P 1/180. ) 

STHI=SIM(TH1«P1 /IBO. ) 

CPHl=CnS(PHl^Pl/lBO.) 

SPH I =S I M PH I *P I /I SO* ) 

DO SI K~1 ,NW 
DO 51 L=i,MZS 
MI = K+(L-U’J'MW 
. !Vj=NH'^MnDt:s 
VT(^^I)=.(O.C,0.0) 

VA=(0.0,0.0) 

VC=( D.0,0.0) 

HO 52 n = l,.NHAX 
CI=C(II) 

X1 = 0* 

X2 = 0. 

Yl = DDW*( II--1 ) ^TP+(K-1)«DWL«TP 
Y2=Y1 

ZA=TP*{ ZL-OOv^+U r-D-IL-n^OWL) 

ZS = ZA 

VD=(0, 0,0*0) 

VB=(0.0,0.0) 

00 53 1 = 1, \‘P 
Z I=-0. 5*AL*^HL*I 
Rl=COS(ZI*Pi /At ) 

2 1= n-l ) =5^HL*TP+ (L“l )*TP=<^ZL/NZ 
Z2=Z1+HL«TP 

YA=( l-l) {K-D ^JP^ZL/Hl 

yr=ya+hl^tp 

CALL 2FFD(X1,Y1,Z1,X2,Y2,22,HL*TP,CTHI ,STHI,CPHI ,SPHI ,VTA ,VPA) 
call ZFPDIXI ,YA ,ZA,X2,YB,Ze,HL*TP,CTH! ,STHI ,CPHI tSPHI ,VTB ,VPB ) 
VB=VB+VTA*EI 
Vn=VD+VT&^Fl 
53 CONTI^UF 

VC=VC+VD«DOW’*C I 
VA=VA + V3*DDW’«'C 1 
52 CONTINUE 

FT(NI) =VA/3./DWL 
FT(NJ)=VC/3*/DWL 

VT(!Mn=ET(Nn /CMPLX(C.O,- 60 .*PI ) 

VT( NJ)=ET(NJ)/CMPLX (0.0,“60*^Pn 
51 CONTINUE 

IF (PHI .FQ. 0.0) CALL CROUT( Z,VT, MODE!, IOM,0, 1,1) 

IF (PHI *FQ. 45.) CALL CRQU T( Z , VT , MOOET^ 1 DM, 0 , 1 , 2) 

EST=(0. 0,0.0) 
on 123 1=1,M0DET 
FST = FST+FT( I ) *VT( 1) 
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123 cr>N’Tl\'Ur 
F.STA=CABS(EST) 

eSTP=rATAN2( AIMAG( EST ) »RFAL( BST) )’«=180 „/Pl 
SIG = 4,*P1’PFSTA^FSTA 
WRITet6,124) ZL?SIG,ESTAjFSTP 

124 FORMAT ( "■»X»4Fi0o4) 

PH I =PH 1 +43 

IF (PHI a>LF«. 43*.0) 0,0 TO 100 
ZKL=ZKX+0»4 
C ND 

COMPL»^X FUMCTION ZMM ( OL »HL , SL ) 

RHAL L^LEfLL 

B=6.2R3ia53 

0=OL 

L=HL 

LF = HL 

HC = SL 

RLF=S*LF 

H=ABS(HC )*L 

LL=LF 

HPL=H+LL 

HP2 L=H+ 2 eO^LL 

HP3L=H+3*0^LL 

HML-H-LL 

SBL=SIM(BLE) 

r.BL=C0S(6LF) 

SBH=SIN(B*HI 
C9H=^C0S( 6«H ) 

SBHML = S 1M(8«HMU 
CBHML = C0S(B’«'HML) 

SBHPL=S IMIB^HPU 
CPHPL=COS(B«HPL I 
SBHP2L=SIN(e*HP2L ) 

CbHP2L = C0S{ B«HP2U 
SBHP3L=SIN( B^HP3L) 

CBHP3L=C0S( B4HP31 ) 

T EMP = SGRT ( D«0+H’^H ) +H 
vi=rad*d/tfwp 

U1=64TFMP 

TFMP=S0RT{D*D+HML^HML)+HML 

UO=B*TEMP 

VO^B^^D^D/TEMP 

TEMP = SORT(n^D+HPL’»'HPL> + HPL 

'J3=8*TFMP 

V3=B«D*0/TFMP 

T FMP=5QRT ( 04^D+HP2L^HP2L )+HP2L 

U2=^B’«'D’4'D/TFMP 

V2=B^TEMP 

TEMP = SOR T(D«n+HP3L’5'HP3U+HP3L 

u4.=:B5SfD*D/tEMP 

V4=iB*TFMP 

CALL SICl (SIU0,CIU0»U0) 

CALL Srci ( SIUI ,C lUl fUl ) 

CALL SIC! ( SI V2,CIV29V2? 
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CALL Siri (SIV4tCIV4,V4) 

CALL SICl (SIU3,CIU3,U3) ‘ 

IF (D.LC^O.O) GO TO 80 
CALL SICI (SlVltClVltVi) 

CALL SICI (S1V0,CIV0,V0) 
call SICI f S IV3,C IV3tV3) 
call SICI (SIU?,Cni2,U2) 
call SICI (S1U4,CIU4,U4) 

.0«(C ?H.ML ^(CIUO+ClVO-ClUi-CIVl )~S8HML^(~S lUO+SI VO +S lUl-S IVI )+C 
2BHPL^(2.«C I V3-^2.’«'C IU3-CIU2-C lV2-CIUl-CIVl)+SBHPL’!'<~SIV3-»-SlU3 + SlU2~ 
35 IV2-S iUl + SI V14-S J U3-$ IV3 ) +CBHP3L’> ( ~C IU2-C I V2+C 1 U4+C I V4 > ■►SE^HP3L* ( S I 
4lf2-S 1 V?-SIU4+$1 V4 ) ’►?*^CRL*CBH«(-CIVl-CiUl + CrV3+C IU3) +2»«CBL=»SBH*< 
3S IV1“S IUl-5 IV3 + S I U3 > -^2. aCBl^CBHP 2L*( C T V3 + C ILf3-C I U2-C I V2 I '►2 **CPL=A $8 
6HP2L’J‘(-SIV3 + S 1U3 + SIU2-S1 V?) ) 

X=rl5,0* fCPHML^(-SlUO-SIVO‘^S^Ul + S 1 VII ~S6HML’#=(“C lUO+C I VO-fC Ufl-C I VI 1 + 
2CPHPL=9^(-2.«SI V3-2.«SIU3 + SIU2 + SIV2 + SIU1+SIV1 I + SBHPL* ( ClV3-»-2.*C 

3 I U3 4-C I U2-C I V2-C 1U1-**C I VI ) 4CBH P 5L^< S IU2 + SI V2-S IU4-S I V4 I ■♦'SBHP 3 (C lU 2 
4-Cl V2-CIU4+CI V4)+2.^C8L^CBH# (S 1V1+SIU1-SIV3-S1U3)+ 2 . ^CBL*SBH* { C 1 V 
51-C rJl-C I V3+C IU3l+2*=^CBL*CeHP2LX^(-SI V3-SILf3-*-SIU2 + SI V2 I + 2 ♦ *C B L^SBHP 
6?L<‘(*-CI V3 'kCIU3 + CIU2-C1V2) ) 

GO TO 90 

ao contimuf 

R = 15.0« fCBHML*(CIU0~r U/l + ALOGIH /HWL I ) +SRHML^ ( S lUO-S lUU + S8HPL* 

2 (2 *^S IU3 -S I V2-S lUl ) +CBHPL’«‘ ( 2 •^C1U3~C 1 V2-C lUl + A LOG (HP2L/HPL ) 
3ALQG(H /HPL) ) '♦'CBHP3LX' { C I V4-C 1 V2 + ALOG ( HP2L /HP3L I I + SRHP5L ♦ ( S I V4-S I V2 

4 ) + 2*'->CBL’5=CBH^ (C 1U3-CUJI+AL0G(H /HPL ) ) +2 . *CBL*$BH’^( S IU3-S I U1 ) + 

3 2* «CBL^C^HP2L* (C 1U3-CI V2 + ALOG(HP2L/HPL) I '♦'2 .*CRLX^SBMP2L«( S IU3-- 
6 SIV2)) 

X=1 3*0* (CBHN!L*( SIUl-S lUO) -►SBHML^IC lUO-C IU1'^AL0G(H«L/H ) ) + CBHPL<' 

2 IS IV2 + SUU~2.«S1U3) + $BHPL 2 .X'C IU3-C I V2-C IU1+AL0G( HP L/HP2L ) 

3AL0G( HP L/H ) I ‘►C3HP3L^{ SI V2-S IV4)4SRHP3L*< C IV4-C I VZ-*- A LOG (HP3L/HP2 L ) 
'►2**CBL*CPH«(SIIJI-SIU3U2,*CBL’*^SBH^(CIU3-C1U1-^AL0G (HPL/H ) 
3L’ACRHP?L«(SlV2~-SIU3)-^?.^CBL«SPHP2L*(ClU3-CIV2+ALaG(HPL/NP2L) I I 
90 ZM'^ = CMPLX(RtX)/(S3LASBL) 

RPfURM 

END 

SUBROUTINF S I C 1 ( S I , C I , X I 
Z=^ARS(X) 

1F<Z-4,0)1,1,4 

1 y=(4.0-2 )#(4.0-*-Z) 

s 1 =^- 1 . 570707^0 

1FIZ)3,?,3 

2 CI=-1.0F-38 
RETURN 

3 S ( < ( ( ( i *7 53 I41F-9 ^y + 1 • 568988E-7) *Y + 1 .37416SE-5 I *7+6 *9 ^9BP9E~4 I 
2=«‘Y+1.964 88 2F-2 )«Y+4.395 509F-1 + SI/X) 

CJ=( <5#772156E— 1>AL0G(Z) )/Z— Z*( ( ( ( (1 •3869R5E— lO^Y**-! • 5 8499 6 E— 8 I*Y 
?+l • 72375 2F-6 ) *Y+ 1 • 1 85999E-4) =^Y+4 .990920F-3 ) *Y+ 1 . 3 15 3 08E-1 1) Z 
RETURN 
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S1 = SIN(2) 

Y=COS(Z) 

Z=4*0/Z 

L; = ( (((((( (A*a48069E~3’J^2“2»2T9143e-2)=frZ + 5a515070F-'?J *2“7.?6l64?F-7) 
2+Z+4,98 7716F“2)*Z-3c33231QE-3)=«‘Z-2,314 617F~2HZ-'1,134<^5PF-4)^2 
3^6.2300 in— 2 ) *Z + 2, 5 83989 £-10 
V- ( ({(((( ( (-3, 108699F“'3*Z + 2*^819179£-2l*Z~6o537283E“2)*2 
? + 7*902034E-2 ) ♦Z“4*400416E~2 K=«^Z-7 ,945 53 6E~3 ) *Z + 2 • 601 293E->2 HZ 
3-3,764 000 6 - 4 ) 12241 8 F~2 ) 64644 lE-7)*Z-^2 o 5 0000 OE-1 

C 1-Z^(S I#V-Y^U) 

5I = -Z’«'{5i]^U + Y*V) 

IF( X) 5j6,6 

5 5 I=r-3, 141593E0-SI 

6 RFTU??N 
Em 

SUBROUliNF CPOUKC, S,N, ID?^, I SYWylWR? 112) 

COMPLEX CnDM9]DM)»SUDM) 

complex F^PjSS^T 

2 FORMAT ( IXv U59lE10,3 tlF15,7^ 1F10,0) 

5 FORMAldHO) 

\ E{ 1 12,^*F,1)GG TO 22 
IF(N,FO.I)S( 1)=S{ 1)/C(1,1) 

IF(N,F0.1)GO TO IGO 
1 F{ I SYM ,N£,0 ) GO TO 8 
00 6 I-UM 
DO 6 J=ItN 
C ( J ? 1 )=C ( 1 , J) 

6 COMTiNUF 
8 CONTINUE 

F=C(l,i) 

DO 10 L-2,N 

10 t (1,L)=C< 1,L)/F 
DO 20 L=2,N 
LLt^L-1 

DO 20 I-LfN 
F=C(I,L) 

DO 11 K:=1,LL£ 

11 F=F-C(I,K)^C{K,U 
cn,u=F 

IF( L.FO, nco TO 20 
P=C(LtL) 

I F( ISYM.EQ.O )G0 TO 15 

F=C(L»I) 

on 12 k=i,lll 

12 F~F-C(L,K)=^C(K,n 
CrL,I)=F/P 

GO TO 20 
15 F=C ( ItL ) 

C(L,n=F/P 
20 CONTINUE 
2 2 CONTINUF 
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DO 30 L=i#N 

P<(L,U 

T:^S{L> 

IF(L.E0.1)GO id 30 
LLL=L-1 
DD 25 KjrltLLL 
?5 T=T-C(L,K)«S(K) 

30 S(U=rT/p 

DO 38 L=^2tN 

1=M_L+1 

11=1^1 

T=S{1) 

DO 3 5 K = II ,N 
35 T=T-C(IrK)vS(K) 

38 S(I)=T 

I F( IWR.LE .0) GO TO 100 
CNOR^.O 
DO AO I=ltN 
SA=CA8S(S(n> 

1 F( SA♦GT.C^^OR ),CNOR=SA 
AO CONTINUE 

I FTCNOR.LE.O. )CNDR=1* 

DO AA 1 = 1 
5?S=SU ) 

SA=CA6S(SS) 

SNDR'SA/C^OR 

p 0 

IF( SA.GT.O.) PH=57.29 578^ATAN2 (AIMAG( SS > ,REAL(SS) ) 

WRl TF ( 6, 2 ) I, SNORtS A, PH 
AA CONTINUE 
WR1TE(6,5) 

100 CONTINUE 
rpturn 

FND 

COMPLEX FUNCTION Z I J ( Xl t X2 t X3 , Y1 »Y2* Y3 t Z1 f Z2tZ3» XA, XB tXCtYAf YB tYC , 
2ZA,ZB,ZC) 

COMPLEX P11,P12,P21,P22 
COMPLEX Qnt012*Q2lfQ22 
COMPLEX SlltSX2f S21tS22 
COMPLEX Pll,R12tR21»R22 
AK=0.005*2*0*3#1A1592 
1NT=0 

012 = S0RTnXl-X2 )*(Xl-X?) + (Yl-Y?)«(Yl-Y?) + {Zl--Z?)’^(2l-"22n 
023 = SQRTUX2~X3)*{X2-X3) + (Y2-Y3)»(Y2-Y3)-MZ2~23)*(Z2-*Z3) ) 

DAB = SORT ( (XA-XB)*( XA-*XE)^“( YA-Y8)«(YA-YB) + ( ZA-ZB)* (ZA-ZB) ) 

DC5=SQR T{ (XC-X6)*(XC-XB) + (YC-'YB)*tYC-YR)<^<2C-'2B>^ (ZC--Z6) ) 
CD1?=C0S(01?) 

SD12=SIN(D12) 

CD23=C0S(023) 

SD23-S1N(D23) 

COAB=COS(DAB) 

50AB=SIN(DA6) 
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CDCB=CaS(DC«) 

S0CB=S1N(DCB) 

CALL ZGS(X1 ,Y1 ^Z1 9X2yY2f XA^YAjZAf XB 9YB9ZB«AK?ni2 9Cni2, S012, 
?DAB ,SDABt !NT»P11 ,P129P21 ,P22) 

call ZGS( XI ,YUZl9X2,Y2,22,XB,YB9ZB^XC9YCpZC,^K,P12i,CD129 SOI?, 
2DCB,SDCB, 1NT?Q11,Q12,Q21 5.022) 

call ZGS ( X2 j Y 2% Z2 V X3 9 Y3 9 Z 3 t X a f Y a , Z a , XB , Y 6 9 ZP p AKY D23 » C D?^ » SD?3 , 
2DaS9SDAB,P^T9Pll,R12,»^2l9R22) 

CALL ZGS( X2 9Y?9 Z2 pX3 pY3 ? Z 3 » XB » Y8 ^ ZB , XC p YC ? ZC 9 AK , D23 p CD23 9 SD23 p 
2DCB pSDCPp iNTYSlltSlZpSZl ,S22) 

Z IJ^P22'^02I+R 12^^511 ’ 

RETURN^ 

FW 

SUBROUTINIF SIMWC(NMAXpC) 

DIYEMSION C(31) 

DO 1 N^IpNMAX 
XMN=FLOAT(N) 

!NIN-N/? 

TT~XNN/2a 

D1F = TT-FL0AT(\»N) 

iMC=2 

IF(DTF^FQ,0<,) NC = A 
I.F( M.FO, 1 ,0R .NaEQ^NHAX) NC = 1 
C(N)=?VC 
1 CONTINUE 
R FTURM 
E m 

SUB ROUT INF Z F FD ( X A 9 YA 9 Z A 9 XB 9 YB * ZB »Df CTHp S TH 9 C PH 9 SPH 9 E T 9 F P ) 
COMPLEX ETpEP ,FS»F JApFjP . 

XAB=XB-XA 
YABsr-YB“YA 
Z A3=^ZR-ZA 
CA = XAR/L> 

CB=YAB/D 

CG-ZAB/D 

CA’>^CPH + CE + SPH)’>STH^CG«‘CTH 

ET=(C.0t0«0) 

FP=(0.090-0) 

1 F( GK ,LU 0.001 ) GO TO 200 
3 =XB’^STH*CPH+Y8^STHYSPH + ZB*CTH 
E JR=CMPLX(CnS(B) ,SINJB) ) 

SKO^STMrn) 

CKD=CQS(0) 

CGD = CDS CG^^'D) 

E S=6 0. ( ,0, 1 .) ^EJB«(CK0~CG0) /SKD/GK 
T= ( C A^CPH+CB’^S PH ) =«^CTH-CG^ STM 
P=-CA*^SPH+CB + CPH 
FT=FS^T 
EP=FS^P 
200 CGNTINUE 
RETURN 
END 
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SmROUrim ZGS(XAtYA,ZA,XB,YB,ZB,Xl,Yl tZl ^XZtYZ t Z2, AK , 
20S,CDS,S05*DTf SDT, INTtPll ,Pl2.P?ltP22) 

CUMPL*^X CST,F J1 tB J2 JA, FJBffFRlf ER2, ETl,ET2,Pllf P12tP21 •P?2fGAM 
COMPLFX SGDSfSGDT 

DATA ‘^TAtOAM, PI/376. 727 t< t3,141*>9/ 

CA=(X2-X1)/0T 
CB=(Y2-*Y1)/DT 
CG=( Z2~Z1)/DT 
CAS=(XB-XA)/DS 
CBS=(YB~YA)/D£ 

CGS-(ZB-ZA) /DS 

CC^C^^CAS+CB’i^rBS + CG^CGS 

I F( ABS<CC).GT. 0 .«>P 7 )GQ TO 200 

SZ= I Xl-XA)+CAS+(Yl*-YA)*CBS+(n“ZA)^CGS 

IF{ INT.FQ.OGO TO 300 

CGDS=COS 

SGOS=CMPLX( .OtSOS) 

SGOT=CMPLX( .0,SDT ) 

INS=2*C lNT/2) 

IF( INS.LT.2) 1MS=2 
IP=INS+1 
DELT=DT/INS 
T = .0 

dsz=cc*delt 

PX1=( . 0 ,. 0 ) 

P12=(.0,.0> 

P21=( 

P2?=(.Ot.O) 

AKS=tAK^AK 

SGN--1. 

DO 100 IN=ltIP 

ZZ1=SZ 

ZZ2=SZ-n£ 

XXZ=X1 + T=^CA-XA~SZ*CAS 
V YZ = Y 1 + T S-Y A-SZ’J^C B S 
7ZZ = Z14-T*CG-Z A-SZ^'CGS 
PS=XXZ*’«^2 + YYZ=^*2 + ZZZ**2 
R1=SQRT(RS+ZZ1**2 ) 
e JAsrCMPLXICOS (R1 ) ,-SIN(P 1 n 
PJ1=FJA/R1 
R 2 =S 0 RT (RS+ZZ 2 +^ 2 ) 
r JB=:CMPLX(C0S(P2) ,-SIN(P2) ) 

FJ2 = EJB/R2. 

ERl =£JA«SGDS+ZZ 1*FJ1*CGDS-ZZ2*FJ2 
FR2=-FJB*5GDS-^ZZ2*EJ2^C&DS--ZZ l^FJl 
P AC = .0 

1 F( RS.GT.AKS ) FAC= (CA*XXZ*^CB*YYZ + CG*ZZZ) /RS 
f.Tl = CC* ( E J2-FJltCGDS) +FAC*ER1 
ET2 = CC^U-J1-E J2»CGDSKFAC^eP2 
C=3.+SGN' 

iPfIN.EQ.l ,0R. 1^I.EQ*IP)C = 1. 
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Cl = C«SlN<<Df*-T) 

C7 = C^‘^P'HT) 

PU = Pll+tTl>fCl 
Pl? = Pi 2 ’^^Tl♦C 2 
P 2 ], = P 2 l + rT 2 «C 1 
P 2 ?=P 22 +FT 2 «C 2 
T=T 4 -oeLT 
SZ=SZ+OSZ 
100 SGN=-5GM 

C ST = "'< «0, K ) =«‘ETA^D,ELT,/{ 1*SGDS*SG0T J 

Pn=CST*PIl 

P 12 =CSr«P 12 

P 2 i=CSTHcPi'l 
P22 = CST=^P22 
P PTIJRN 

200 SZ1 = ( Xl-XA ) ^CAS->-( Yl-YA) ^CBS+ ( Z 1-ZA HCGS 

RHl^SORT ( (X l“XA-SZl=«fCAS ) ( Yl- Y A-S Z 1*CP S ) >4^* 2> { 71--Z A*-FZ 1 ) **2) 

SZ2 = SZl-H>T«CC 

PH2 = S0RT( (X2-XA-SZ2+CAS )*=^2>( Y2“YA“SZ?’»‘CBS > Z2-ZA-$Z2*CG^)^^'>) 

DDY = (pHl-»-PH2>/2* 

IFfDDK»LT^AK)DOK=AK 

C ALL 2GKW( *0,DS?SZ1 ,SZ2,DDK,CDS,SDS^$DTyI o ? Pll 9 P 1. 2» P2 U P2 ? ) 


300 SS=SQRT( l^-CC’t'CC) 

C AD= {CGS^CP-C BS^CG) /SS 
C«D - (CAS*CG-CGS>^CA) / SS 
CGD* (CBS*CA-CAS^CB)/SS 

OK- ( Xl-XA) ^CAD'f {Yl-YA)ACBD+( Zl-ZAJ’^CGD 
DK^A8S( OK) 

I P(DK.LT, AK )OK=AK 
XZ=XA4SZ^CAS 
Y2-YA+S2*€BS 
ZZ = ZA>SZ*CGS 


XP1=X1~DK^CAD 

YPl=Yl~DK»^CBr> 

ZP1=ZI~0K«CGD 
C AP=C BS*r.CD-«CGS+CBO 
CPP=CGS*CAD-CAS’i^CGD 
CGP=^CAS^CBD“CBS^CAD 

P l-CAP* ( XPl-XZ)+CeP*( YPl-YZ) +CGP*< ZPl-ZZ ) 
T1=P1/SS 




CALL ZGMMiSlrSl+DS^TltTl+nT 

R FT URN 

FND 


9DK9CDS9SOS 9SDT9CC9PI 1 9P 12 9P2I tP22T 


SUBPOUTINF EXPJ(\/l,V2,Wi2) 

COMPLEX PC »F159 S? TvUCtVCtVl 9 V2, WI 29 Z 
DIMENSION V(21) 9W(2X) 90(16) ,8(16) 
data V/ 0*22284667F 00, 

20.1 18R9321F: 01,0.299273638 C? 1 , 0 .577514E6E 
20. 15 9828748 02,0.933078128-01,0.492691768 
20.226994958 01,0.366762276 01,0.542533668 
20.1012022^^8 02,0.131302628 02,0.166544088 
20.25623894E 02,0.314075198 02 , 0 . 38530683E 


01 ,0.9837 46 74E 
00,0.121559546 

01.0. 756591626 

02.0. 207764799 
02,0o4e0260n66 


01, 

01 9 
01 , 
02 , 

02 / 
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DATA W/ 0.A5896460T 00, 

20 ,41700083^ 00,0. 11337338P 00 ,0. 103991 97F-01 ,0 . 2610 1720^ -03, 
20.P9884791F-06,0.218234B7E 00 ,0. 342210 17E 00 , 0 .26302788 E 00, 
20.1 2642 8 82 F 00 , 0 .40206B65F-0 1 , 0. 85638778E-02 ,0.121243610-02, 
20.1 1 16 74408-03, 0. 645992 67E-03,0. 22263 169E -06, 0.42274304F-08, 
20.5 9218 973E-1 0, C. 14565 162 E-1 2 ,0 .148302708-15,0.160059498-10/ 
DATA D/ 0.22495842E 02, 

2 0.7441 156SE 02,-0.414315768 
2 0.16021761^ 03,-0,238621958 
2 0.122547788 02,-0.101619768 
2-0,210695748 02, 0.220464908 
DATA E/ C,211C3107E 02, 

2-0.3795978TE 03,-0.974892208 
2-0 . 1 29 1 0 93 1 E 03 , -0 . 5 57 0 5 574 E 
2 0.17949528P 02 ,-0 . 3298 10 14E 
2 0.22236961E 02, 0.39124892F 
Z=V1 

DO 100 JIM=1,2 
X=^9EAt(2) 

YsAlMAG(Z) 

E15=(.0,.0) 

AB=CABS< Z ) 

I F( AP.EO.O. )G0 TO 90 
IE(X.GE.O. *A\'D. Afi.GT* 10. )G0 TO 80 
YA=ABS(Y) 

IF(X.LE.0. .A?^D. YA.GT.IO.IGQ TO 80 

IF( YA-X.GE.17.5.0R.YA.GE.6.5.0R.X^YA.GE.5.5.0R.X.GE.3.)G0 TO 20 
1 F( X.LE .-9.) GO TO 40 
IF( YA-X.GE.2.5)GO TO 50 
1 F( X + YA .GE.l .5>GQ TO 30 

10 N^6.+3.*AB 

E15 = l./<N-l. )-Z/N«’«^2 

15 N=N-i 

E 15 5^ 1 . / ( N-1 . ) -Z *E 15/N 
1F(N.GE.3)G0 TO 15 

E15=2«E15-CMPLX( .577216+AL0G( AB) ,ATAN2(Y,Xn 
GO to 90 

20 Jl=l 
J 2= 6 

GO TO 31 

30 Jl = 7 
S2-ZI 

31 $:=(.0,.0) 

YS=Y*Y 

no 32 r=ji,j 2 

XJ = VM)+X 

CE=W( I) /(XITXl+YS) 

32 S=S+CMPLX(X1TCF,-YA*CF) 

GO TO 64 

40 T3=X*X-Y*Y 
T4=2.*X’«=YA 
T5=X=^T3-YA=«fT4 
T6=X*T4+VA»T3 


03, 

-0.787543398 

02, 

0.11254744E 

02, 

03, 

-0.50094687E 

03, 

-0.6648 78 54E 

02, 

02, 

-0.47219591E 

01, 

0.79729681E 

01 , 

01, 

0.89728244E 

01/ 



02, 

0.12900672E 

03, 

C. 179492 26E 

02, 

03 , 

0.13524801E 

02 , 

0.14696721E 

03, 

00, 

0.31028836F 

02, 

0.81 6576 57E 

01, 

02, 

0.81636799F 

01/ 
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UC=CMPLX<D< 11 ) 13 ) *T3«J-T5--E « 1 2 ) ’J'YA-F (13 ) 

2 ^'(H)-^F(12}^X-^F(13)+T3-6^T6 + D(1-2J^YA^D(13)^T4I 

VC = CMPLX ( D< 14 ) ( 15 ^ *X + D(16)^T3^T5“E C 1 5 ) * Y A~F U 6 ) «T4 , 

2 F ( 14} '►F ( 15)*X4^F( 16) >^' T3■^T6-^D ( I 5 ) ’f^YA'^D ( 16 )*T4) 

GO TO 52 
50 T3=XXfX-v=«>v 
T4=2 ,^X>f^YA 
T5=X«T3-YA«T4 


5 2 


54 


56 


80 


PO 

100 


76=X*74+YA^T3 
T7=X»5^T5~YA«T6 
lB = X’i‘T64YA#7 5 
TP=X^T7-yA«Te 
TlOrX + Te-i^YAAT? 

UC = CMPLX( DU ) *t-D( 2 3) ♦T3-«-D( 4)^754-0 t5}^T7^19-i F ( 2 ) *Yfl'^ F ( ^ 

2 .F 4 *T6*P(-i)»T8),fm*E(2)*X*E(3)*T3*E(4)*T5 + Ef5)lT7+T10r 

3( D( 2)*YA + D(3> ♦74+DJ 4)’>T6 + DS 5 )’5'T8 ) ) 

VC~CMPLX(D(6)+D(7) *X + 0( 8) AT3+DC 9)*75+D ( 10 ) ^T7•^T9-( F ff 7)<tYA^F ( P 74 

?,F 9 *T6*ni0)*TB),6(6)^t(7)*X.E(8)4T3^E(9liT5iF(10)*T7lnCK 

,(D(7)*YA + D( 8 ) *T4+D (9)’^T6-»^D(10UT8)) 

FC=^UC/VC 

S=FC/CMPLX( X 9 YA ) 

FX=EXP(-X) 

T=:EX*CMPLX(CO.F (YA) ,-SlN(YAU 
FI5=S«T 

IF(Y.L7i,0, ) F15 = CDMJG ( E15 ) 

GO TO PO 

2+2. <^6788)^ 

^ (7+8^18215) + ol5P654F~4 / (7 

U?.7342)<-. 3170318-7/(7+19, 3957) ^/(/ + 

E15:=F15^CFXP|-2 ) 

I F( J IM, FQ, 1)W12=F15 

7 —V 2 

7=^V2/v1 

TM=aTAM 2( AIMAG( Z) sRFALt 2 ) ) -AT AN? ( A IM AG ( V 2 ) ?RFAL( V?) ) 

2 + ATAN2 ( A 1HAG( VI) «RF AUVl ! ) 

AE=ABS(7H> 

IF( AB.L7»1 ^) TH=,0 
lF(7HeGT,l,)TH-6*2e31853 
1 F ( 7H*LT« — 1 a ) TH = -6&2B31 853 
W12 = W12-F15+Ci^PLX( .OUH) 

RETURN 


F NO 

LUNPLFX Ff^, ,2)fF(^,2),GAM,Pll^P12yP21,P22 

COTPLfJ LNfO";D!l:!!^F'!'’?»’ESl,ESP,FTl,F,T?,FXPA,FXPB 

COMPLEX 


fb»ec^fk 

FGZ(2t2 


eXA( 2) ♦FXB(2) 
data fta,gam 
DSQ*D:«^n 


'-»£»cMisAn,Kii^Pl^,P21fP22 

KfFL tFKL , FG2I , ESI, 852^8 71, FT2 ,FXPA 

) ?GM (2) ,GP (2 ) 


?P2?) 


P(2) 

»PI/376*727, { .0,1 


1 «3.14l«;6/ 


SGDS=SGQ1 

IF{ SZ.LT.Sl )SGDS--SGD1 
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5G0T=SG02 

IF( T2.LT.T1 )SGDT=«SGD2 

IF( ABS(CPSl) ♦GT.0.997)G0 TO 110 

FS1=CEXP(GAM«S1 ) 

FS2=CFXP (GAM*S2) 

FT1=CFXP{GAM^T] ) 

FT2=CEXP {GAM#T2) 

C^O/SOftTU .-CPS I *C PS I ) 

B=C^CPSI 

ER=CEXP(GAM*CMPLXC .OtS ) ) 

FC=CFXP (GAp}«CMPLXC.O,C) ) 

DO 10 K=l*2 
DO 10 L=l»2 
10 F(KtU = ( .0,.0) 

EK = EB 

DO 50 K=l,2 

FK=(-1)+^K 

FL=FC 

DO 40 L=l,2 
FL~(-U^*L 
eKL=EK^EL 
XX=FK*B+FL*C 
S1 = S1 

DO 30 1=^1, 2 

P1~S0RT(DS0+SI’^SI+T1*TI-2**S1*T1*CPSI ) 
R2=SqRT(DS0+SI*SI+T2*T2-2.*Sl*T2*CPSI) 

CALL FXPJ(GAM*CMPLX(R l+FK»Sl+FL*Tlt-XX) , 

2 GAM«CMPLX(R2+FK*SI+FL*T2,-XX) tEXACI ) ) 

CALL EXPJ(GAM’^CMPLX(R1+FK^SI + FL’0^T1,XX) , 

2 GAM*CMPLX(R2-fFK*Sl-i-FLA^T2fXX) tFXBC n ) 

lFrK.FO.2 .OR. L.EQ*2)G0 TO 30 
ZC==SI*CPSI 
EGZ1=CEXP(GAM^2C) 

CALL EXPJ(GAM«(RI>ZC-T1) »GAM*(R2+ZC--T?) tFXPB) 

CALL EXPJ(GAM=i^(Rl-ZC + Tl) tGAM’MRZ-ZC'fTZ) ,EXPA) 

F (I » 1 )=2.^S6DS*( .0,1. )*EXPA/EGZI 
F(I »2)=2**SGDS*( *0,1.)*EXPB*EGZ1 
30 SI=S2 

E(K,L)=E(K,L)+(EXA(2)~EXA(1) )*EKL^(EXBC2)-EXB(1) )/FKL 
40 EL=1./EC 
50 EK=1./EB 

CST=-FTA/< 16.^PI*SGDS*SG0T) 

P11=CST^{ ( F( 1, 1)+EC2,2)*ES2~E(1 t2)/ES2|*ET2 
A +(~F(lt2)-E(2,l)*ES2+E(l»l)/ES2)/FT2) 

P12=CST*( (-F( If U-E( 2,2) *ES2-»-e( 1,2) /E$2)«ET1 
B ♦{ F(l,2)+EC2,l)«eS2-E(l,l)/FS2)/ETl) 

P21=CST=*( (-F( 2, 1)-E(2,2) *ES1^E(1,2)/ES1)*ET? 

C F(2,2)+E(2,l )#ES1-E(1,1)/ES1)/ET2) 

P22=CST^( ( F{ 2 ,n+E( 2 , 2 )^ES 1 -E(l, 2 )/ES 1 )^ET 1 
D ■*'t^F(2,2)-E(2,l)*ESl4E(l,l)/ESl)/ETl) 

RETURN 
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no fFtCHSl •LToO.)GO TO 120 

TA=T1 
1 T 2 
GO T(? 130 
120 TA=~n 
TB=-T2 
SOOT=~SGDT 
130 SI=S1 

no X30 I:=l,2 
r J-TA 

on j-r ,2 

Z 1J=7 J“S1 

P>=SQ^T(r)S0-»-71 J*Z1J) 

W = ^ 4 21J 

f FC2IJ.LT,0*)W^nSQ/(R~2 IJ) 

V=R-7IJ 

IFf ZIJ.GT.Oo ) V=DSQ/ (R+Z IJ ) 

I F ( J , r 0 , 1 ) V 1 

1 FT M'Jl-W 

FOZ< It J)-CFXP(GAH>«'7I J) 

140 TJ=T& 

CALL FXPJlGAr^i^Vl tGA^^^V,GP( 1 ) ) 

CALL FXPJiGAM’^Wl tGAM*W$GM( 1 ) ) 

150 S1=S2 

C ST = f- TA/( e .^P I«SGOS*SGDT ) 
c^ll-CST’i^(GM(2)’»=FG7 ^ 2t2)+GP( 2) /FGZ< 2^2) 
2-CG0C*{GM? 1 )*'«'eG2( I 9 2 ) ^-GP U ) / F GZ C 1 9 2 ) ) ) 
Pl?. = CST^l-GM{2)*FGZ(2tl )~GP ^ 2 ) /FGZ « 2 a ) 
24fX.DS’^ ( G!wi< 1 )*FGZ ( 1 ) -^-G P ( 1 > /FGZ Ut 1 H ) 

«^2l=CST=<'(GMn )>>FGZ( 1,2)+GP( 1) /FGZU,?I 
?-CGn.S’MGM<2 I’FFGZ C ?9?)^-GPr 21/FGZ(29 2) ) ) 
P2?=CbT U^EGZ Cl a ) -GPCl ) /LGZ (1 a ) 

?*»‘CG0S^7GM(?)«FGZ( 2a )^GP (2)/FGZ(29l) } ) 
RCTURN 
F W 
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